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Abstract: We consider here the problem of extrema for the Kac functional with long range, 

^ \ ferromagnetic interaction, and vorticity conditions at infinity which make it not weakly closed. 

! Using a gradient- flow dynamics, we investigate local minima, showing strong analogies with 

^ ■ the Ginzburg-Landau functional in infinite volume. 
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^ 0. INTRODUCTION. 

S: 

I— We consider here the XY or "planar rotator" model and its corresponding continuous 

■ version, with internal continuous symmetry group 0~^{2) and long range, ferromagnetic inter- 
\Q [ action. Such interactions were introduced by Lebowitz and Penrose in Statistical Physics as 
^ [ a generalization of the celebrated Kac, Uhlenbeck and Hemmer model, accounting for liquid 
^ . vapour phase transitions, and giving mathematical foundations to Van der Waals theory. 
O ■ A theorem of Dobrushin & Shlosman [Si, p. 78] asserts that at any inverse temperature 
^ I (3, every infinite volume Gibbs state corresponding to a fairly general, translation invariant 

', hamiltonian H with internal continuous symmetry group 0~^{2) on the 2-d lattice, is invariant 

^ ■ (up to conjugation of charge) under the group 0~^{2). This is absence of breakdown of 

■ continuous symmetry. Later Bricmont, Fontaine and Landau discovered that this Gibbs 
state is unique. It presents many interesting, and yet not fully understood features, such as a 
particular form of phase transition at low temperature, which is characterized by the change 
of behavior in the correlation functions. For the XY system they were described by Kosterlitz 
& Thouless in term of topological vortices. 

Here we look instead at the free energy functional, which provides, through the mean 
field approximation, a good approximation to the canonical Gibbs measure on the set of 
magnetizations m, i.e. suitables averages of spins over mesoscopic regions. Moreover, vortices 
can be enhanced by imposing some boundary condition on H at infinity. 

One usually consider the vector spin hamiltonian on as a limit of hamiltonians (the 
thermodynamical limit) on finite lattices A T? , with boundary conditions on A^. Then 
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the free energy functional takes the form 

T{m\m'') = - J dr J dr' J{r — r')\m{r) — m{r 



(0.1) 



+ - f dr f dr' J{r — r')\m{r) — m{r')\'^ + j dr{^fp{m{r)) — fi3{mj3)) 



where ffj{m) = — ||mp + ^/(m) is the free energy for the mean field approximation, and /(m) 
denotes entropy function. Actually formula (0.1) holds also (formally) in the continuous case, 
see Appendix C. Systems of spins valued in {+1, —1}, i.e. the scalar case, with long range 
interaction are well understood, due in particular to a series of papers by Cassandro, DeMasi, 
Presutti and their collaborators, who studied in great detail interfaces and equilibrium shapes. 

In an attempt to generalize this theory to the XY model [El-BoRo], the authors have 
looked for local minima of the free energy functional, and observed in numerical simulations, 
among other things, vortex configurations on finite lattices A, induced by the vorticity on 
the boundary A'^, which are very similar to those arising in solutions of Ginzburg-Landau 
equations, describing a superfluid or a supraconductor subject to a magnetic field. Related 
examples include the 't Hooft-Polyakov monopole and the Skyrme model. 

We will consider in this paper a continuous version of the rotator, that complies to the 
methods of Functional Analysis. The minimization problem for the free energy functional 
(0.1) and its continuous version in infinite volume has the 0^{2) symmetry, which we break by 
imposing on the magnetization m a vorticity condition at infinity of the form m{x) = mjse^^^ , 
n e N \ 0. Here g]0, 1[ is the critical value for the mean field free energy, and we work at 
a sub-critical temperature, i.e. (3 > 2. 

The obvious conjecture about a minimizer is that there is a radially symmetric vortex of 
degree n, expressed in polar coordinates as m{x) = Un{r)e^'^^ , for some nonnegative function 
Unir) with Un{r) mp as r — >^ oo, and Un{0) = 0. For n = 1, this is called a hedgehog. 

A similar conjecture holds in the case of Ginzburg-Landau equation in a disc, where it 
is known that the hedgehog i/ji, uniquely defined, is stable, i.e. all the eigenvalues of the 
self-adjoint second variation operator L^j , 

1 d'^ 

(0.2) {'v\L4,^v) = - -^8 (V'l + ev)\s=o 

are positive. This was obtained independently by Lieb and Loss [LiLo], and Mironescu [M]. 

In infinite volume, the situation is more subtle, due to the translational symmetry of 
the problem. In fact, Ovchinnikov and Sigal showed (still for Ginzburg-Landau equation) 
that for any n, there is a unique radially symmetric vortex -i/'n of degree n, that minimizes 
^ren('0) among all functions of the form tp{x) = Un{r)e^'^^ . Here £ren is a renormalized energy 
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functional. The self-adjoint second variation operator L^^ is decomposed as channel operators 
in Fourier modes ; the n-th mode is a positive operator, and has not for an eigenvalue, but 
is an eigenvalue for the n + 1-th mode (due to the translational symmetry) , and there are 
also negative eigenvalues for some higher modes. Thus we cannot conclude in this case to 
linear stability of the radially symmetric vortex. 

Minimizers of the free energy of the scalar Kac model, and their stability were investi- 
gated by DeMasi, Presutti et al. In the 1-d case, they proved existence of a solution, called 
instanton, unique modulo translations, of the minimization problem subject to the condition 
that the magnetizations m{x) tend to ±m^ as x tend to infinity. Their result was rederived 
lateron by Alberti and Bellettini [AlBe], and extended to higher dimensions, for a class of 
functions "varying only in a direction e". In fact, both functionals (Kac and Ginzburg- 
Landau) are not convex, and the free energy of Kac model is not even local, so that direct 
methods [Da] don't apply here. 

Rearrangements methods were used in both [LiLo] for Ginzburg-Landau equation in 
a disc, and [AlBe] for the scalar Kac model. Actually, only partial convexity is achieved 
when restricting to the class of radially symmetric functions for Ginzburg-Landau, while 
for the scalar Kac model, convexity is retrieved, by rearrangements, on a set of increasing 
functions. But it is hard to figure out at least what rearrangements would mean in the 
vector (or complex) spin model. We follow here an alternative route, elaborated in the papers 
[AlBeCasPr] , [DeM] , [DeMOrPrTr] , and culminating in [Pr] , which consists in looking instead 
for extrema, verifying Euler-Lagrange equation, as limiting orbits for a certain dynamics, 
known in that context as the gradient flow dynamics Tt (see (1.11) ). Free energy is a 
Lyapunov function for the gradient flow dynamics, and thus we can resort on methods used 
in parabolic equations. 

Our main results are related to existence of radially symmetric solutions nin (x) for Euler- 
Lagrange equations Fi{m) = (see (1.8) ) associated to a suitable renormalized free energy 
functional in infinite volume J-^en defined in Proposition b.5. We have : 

Theorem 0.1: Assume that J is on negative definite in the sense of [FrTo], and the convolu- 
tion operator expressed in the (e„)-sector enjoys asymptotic properties described in Definition 
1.1. Let m{x) = {mp + 'f;(r))e*^^, with v e (see (1.28) ) such that J^ren{m) < oo. Then 
any limit point m* of T^m, t — > oo (in the sense of uniform convergence on the compact sets) 

satisfies Fi(m*) = and .Fren("7.*) < .Fren("7.). 

Thus we can describe some local minima in the space of all configurations m. Next we 
study the spectrum of the n-th Fourier mode of the second variation operator L^* around 
m*. Here we assume a particular form for J. 

Theorem 0.2: Assume J is a Gaussian, and moreover that the critical point m* verifies 
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m*{x) = + v(r))e*"^, with v* e X^. Then the n-th Fourier mode of the second variation 
operator L^* (due to breaking the gauge invariance) has purely continuous spectrum. 

(The decay property on v* could certainly be removed, see Sect. 3. ) So we recover the 
result of [OvSi] for the Ginzburg-Landau functional relative to the n-th mode only. Contin- 
uous,, but also negative spectrum of Lm* near suggests that only "linear instability" (not 
exponential) may occur, in this mode, around the equilibrium. Nevertheless our conclusion 
remains much weaker than for the Ginzburg-Landau functional [OvSi] , since we have no infor- 
mation about higher modes. Naive intuition suggests that is an eigenvalue for the n + 1-th 
Fourier mode, corresponding to the breaking of the translation group. So the "hedgehog 
conjecture" described above in the case of Ginzburg-Landau equation, remains largely open 
in case of the non local Kac functional (here the case n = 1 doesn't seem to play any special 
role) . 

1. The gradient-flow dynamics. 

We look for a solution of the free energy variational problem, among all configurations 
with given degree n > 1, obtained as a limiting orbit of the gradient- flow dynamics. The 
minimization problem in infinite volume involves an approximation process in finite boxes A, 
and renormalization of the free energy. Definitions of the thermodynamic functions are given 
in Appendix C, and we refer to [El-BoRo] for details. 

a) The infinite volume gradient flow dynamics. 

We start to construct the gradient flow dynamics in subject to a vorticity condition 
at infinity. It is convenient to express everything in Fourier modes, using the identification 
L2(R2) Ri ®j^zejL^(R+;rdr), ej{9) = e'^^ . In Appendix B we discuss some properties of 
the degree, and show that higher harmonics can be involved in the Fourier expansion of m, 
provided they decay sufficiently fast as r — > oo ; but because of non linearity, we choose a 
single component j = n, the (cn) sector. 

For given n, Euler-Lagrange equation for J-ren given by (b.28) is the same as we would 
obtain formally, i.e. without renormalizing the energy. It follows that it is of the form 
Fi{m) = 0, where 

1 ^ TTl 

(1.8) Fi{m) ^ -J *m+ -I'{\m\)j^ 

(recall the notation /(|to|) 
to Fo{m) — 0, where 

(1.9) 



/(m) whenever / is rotation invariant. ) Eqn.(1.8) is equivalent 



J * m 

Fo(m) — —m + f{P\J * m\)- 



J * m 
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Here / = = Ii/Iq, and Ii, denotes the modified Bessel function of order v. Solutions 

of (1.8) or (1.9) tend to cluster near the manifold \m\ = mp, which is the critical set for the 
free energy of mean field. Here > for /? > 2 satisfies mp = f{Pmp). 

Let also .Ftt(^) = Jj^^a e~^^^u(x)dx and Hnu(p) = Jq°° Jn(rp)u(r)rdr denote Fourier and 
Hankel transformation, respectively. Hn is defined on the core Cq°(R"'") of smooth, compactly 
supported functions on the half line. For u e Co°(R"''), we have Hnu{p) = 0{p'^) as p — > 0, 
and Hnu{p) = 0{p~°^) as p — > oo. It extends [Ti] to a unitary operator on L^(R;rcir), and 

The convolution operator thus becomes J* = [HnJ^JHn)^^^- We look for a radially 
symmetric solution of (1.9), of the form m{x) = tt(r)e*"'^, with u{r) — > mp as r — > oo, and 
tt(0) = 0. So we get, formally, —AnU + ^/'(tt) = 0, An = HnJ^JHn, or 

(1.10) -u + f{(3AnU) = 

Unless /? < 2, in which case tt = is the unique solution of (1.10), and m = the unique 
minimizer of J^ren (see [El-BoRo,Proposition 3.5] for the case of a lattice, ) the map u i— > 
f{f3Anu) is not a contraction, so following [Pr], we are led to consider the associated gradient- 
flow dynamics, and solve the "heat equation" 

(1.11) = ~" + /(/5^n«), w(0, r) = Uo 

Although An is a bounded operator on L^(R+;r(ir), we are faced with the problem that 
u ^ L^(R+; rdr) because of the condition u{r) — > m/3 as r — > 00. We need the property, that 
An acts naturally upon functions having asymptotics near infinity. It should also map the 
constant function 1 to itself modulo Lp'. So we introduce the : 

Definition 1.1: We say that A^ has the asymptotic property ifi', given x a smooth cut-off 
equal to 1 near infinity, for all k G ^Z, An{-)~^x is a smooth function on the half-line, with 
A^((-)~^x)(r) ~ T~^{l-\-^ + ^ + - • ■), in the sense of asymptotic sums, as r ^ 00. Moreover, 
if Oi = we say that An has the asymptotic property with vanishing subprincipal symbol. 

We don't know if Definition 1.1 could actually be derived from similar arguments in 
Hankel's theorem leading to the inversion of Hn (see [Wa,p.458]. ) Operator An with the 
asymptotic property reminds us of a Pseudo-Difi^erential Operator, with r as the fiber variable, 
and "hidden" phase coordinates that show up when writing the integral representation of 
Bessel functions. It is non-local (not only because convolution is not local, but more seriously 
because a change to polar coordinates always destroys locality, ) with a priori only weak 
decoupling from to 00, see Appendix A. 

Since / is only defined on R"*", operator An must be also positivity preserving, i.e. 
AnU > almost everywhere (a.e.) if u > 0. 
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Example 1: J is a gaussian normalized in L^, so that J- J{p) = exp[— pp^], p > ; we know 
[Wa,p.395] that An{r,r') is given by Weber second exponential integral 

f°° 1 
(1.12) / exp[-pp'']Jn{rp)Jn{r'p)pdp = — exp[-{r'' + r'^)/4p]In{rr' /2p) > 

Jo 2p 

with equality only if rr' = 0. So is clearly positivity preserving (more precisely, positivity 
improving, i.e. A^u > a.e. if w > 0, not identically 0, ) and stationary phase arguments, 
together with the asymptotic expansion of modified Bessel function 



(1.13) In{x) = -= (1 - (n^ - l/4)/(2x) + •••), x^oo 

\' Zttx 

show that r'^An{{-)~'^x)if) ~ 1 + (/c^ — v^^jr^ -\ , so has the asymptotic property with 

vanishing subprincipal symbol. 

Example 2: J^J{p) = e~^^{l +pp), and An{r,r') can be computed taking derivatives with 
respect to p, from the expression [Wa,p.389] 

(1-14) / Jn{rp)Jn{r'p)pdp ^ —=Qr,-y2{ Of ) 

Jo P TTvrr zrr 

where Qn-i/2 is a Legendre function of second kind, for which we have the integral represen- 
tation Qn- 1/2 (cosh r/) = J^(2coshs — 2cosh?7)~-'^/^e~"'*cis, see [Le,p.l74]. It is easy to check 
again that An{r,r') > 0. 

Note that these J's are non negative definite in the sense of [FrTo], i.e. both J and 
J^J are > 0, this excludes antiferromagnetic potentials, or potentials with non definite sign. 
Of course, > implies that An is a positive operator in the mean, i.e {Anv\v) > for 
all V e L2(R+,rdr), but this is not sufficient for our purposes. The additional condition 
(J^J)'(O) = seems to ensure that An has the asymptotic property with vanishing subprin- 
cipal symbol. These features are suggested again by formal stationary phase expansions in 
the integral representation of Bessel functions, but from a ill-behaved phase that becomes 
rapidly oscillating for large r. 

Our first result deals with existence and regularity of solutions of (1.10) in or in 
Sobolev spaces. This will not play the most important role in the sequel, but this is a very 
natural property. We introduce the closed convex set W = {v & L^{'R,'^ ,rdr) : v{r) + rrip > 
a.e.} and W{[0,T]) = C^{[0,T],W) with norm 

(1.18) \\vit, r);Wi[0, T])\\ = sup,e[o,T] \Ht, r) 

Note first that, since J > and / J = 1, we have J^J{p) < 1 and ||^n|| < 1; where || • || 
denotes the (C*-) norm of operators. 
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Theorem 1.2: Let as above An be positivity preserving, and enjoy the asymptotic property 
of Definition 1.1 with vanishing subprincipal symbol. Then for all T > and all vq G W, 
there is a unique v e VF([0, T]) such that u = v + verifies (1.11) with tt|t=o = vq + mp. 

Proof. Since A^ is positivity preserving, AnU > a.e. if tt > a.e. so (1.11) makes sense 
foTu = v + mp, V e W. Moreover, A^u = A^v + vi + mp, Vi = m/3(^„l — 1) £ L^(R+) 
because of the asymptotic property with vanishing subprincipal symbol, and of Proposition 
a.l. Equation (1.11) with initial condition v{0,r) = vo{r) can be written in the integrated 
form 

(1.19) v{t, r) = e-'vo{r) + f dtie''-y{pAnv{h, r) + Pvi{r) + (3mp) - f{(3mp)] 

Jo 

Let again w{t, r) = v{t, r) — e~*vo{r), this rewrites as w{t, r) = ^{t, r, w), with 

(1.20) ^t,r,w)= I dtie''-y{pAnw{tur) + Pmp + Pvi{r) + Pe-''AnVo{r))-f{Pmf3)] 

Jo 

Given Vq e W, denote by Wt,Q([0, T]) the closed convex set of C°([0,T],L^) consisting of 
all functions w(t,r) such that w(t,r) + e~^vo(r) e C^([0,T],W). Because / > 0, $ maps 
Wi,q([0,T]) into itself. Moreover, since ||/'||oo = 1/2, the estimate 

(1.21) mt,r,wi)-^t,r,W2)\<'^[ dhe''-*\An{wi - W2){ti,r)\ 

^ Jo 

and the fact that ||^n|| <lonW prove that if /?^T(1 — e~^^) < 8, then $ is a contraction 
on Wi,q([0,T]). Theorem 1.2 then follows from the group property. J^. 

Consider next the derivatives of u. Applying the radial field rd^ to J*m{x) — e^'^^Anu{r), 
X = xi + 1X2 — re^^ we find, in the distributional sense : 

d dJ dJ 

(1.22) e'-'^^r—Anuir) = xi- — * mix) + a;2T; — * mix) 

or axi 0x2 

which shows, since u = \m\ : 

(1.23) \drAnu{t,r)\ < ||VJ||i||w||oo 

II • 111 being the norm. On the other hand, using (c.6) we compute 

(1.24) ^ +-)An = Bn, Bn{r,r')= / dppJn-i{rp)Jn{r'p)TJ{p) 

or r' Jo 

where as before Bn{r, r') stands for the kernel of operator Bn with measure rdr. We have 
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Lemma 1.3: With notations above, < 1, ||^^n|| < for some C > 0. 

Proof. The first inequahty results from B^Bn — Hn{J^J)^Hn and the fact that is unitary 
on L^(R+,r(ir). We use also (a.l4) to write 

Since the multiplication by pJ^J{p) is bounded on we get the second inequality. X 

From this we can extend Theorem 1.2 to show regularity of the r-derivative u' of w. 
Namely, let be the closed convex set {v & W : v{r) + mp < 1 a.e., v' e L^(R"'", rc/r)}, 
and W?^([0,T]) = C°([0, T]; W?^) with Sobolev norm as in (1.18). We have : 

Proposition 1.4: Let v e C°(R"'", W) be the solution of (1.11) constructed in Theorem 1.2, 
with initial value vq^W'^. Then for all T > 0, v e W^{[Q, T]). 

Proof: With notations as in the proof of Theorem 1.2, we have 

— (t,r)=/3 / dse'-'f{l3Anw{s,r) + l3mp+l3vi{r)+l3e-'AnVo{r)) 

(1.25) -^0 

d 

X —{Anw{s,r) + vi{r) + e~''AnVQ{r)) 

We denote by $'(t, r, zw) the RHS of (1.25) and prove that $' is a contraction on VF^([0, T]) if 
T > is small enough. But this results also from Lemma 1.3, and the fact that w e VF([0, T]). 

This Proposition extends easily by induction to all derivatives, so Sobolev embedding 
theorem shows that r ^ v{t,r) inherits the regularity of its initial datum. 

Existence and uniqueness result in however, falls far short of our needs to ensure 
existence of a limiting orbit satisfying (a.lO) as t ^ cxd, or to provide suitable asymptotics, 
essentially because we lack of a uniform bound on \\v{t, r); VF([0, T])|| as T ^ oo. We restrict 
henceforth to continuous functions that tend to sufficiently fast as r —>■ oo. Such initial 
conditions will be used in the Barrier Lemma below. It turns out that the evolution equation 
doesn't either provide a uniform bound on the norm of u{t, r) but we can still obtain 
indirectly such estimates. Asymptotic property of operators An will be used as a hint to 
model our functional spaces. So we consider the space of functions v G Co(R"'"), such 
that sup^g[i |r'^t'(r)| < oo. It is easy to see that is a separable Banach space, with 
norm 

(1-28) ll^^;^fell =sup^e[o,i] 1^^(01 + sup^e[i,+oo[k''f^(OI 
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For A; = 2, we have C L^. In Appendix A, we give continuity properties of An acting 
on these spaces. Let also X° be the closed convex subset of consisting of functions 

V e Co(R+), such that v(r)+m/3 > 0, and define X^([0,T]) = C°([0, T], X^). 

Theorem 1.5: With the assumptions above, and if in addition has vanishing principal 
symbol, then for all T > small enough, and all vq G X^, k = 1,2, there is a unique 

V G X^([0, T]) such that u — v + mp verifies (1-11) with u\t=Q = f o + 

Proof: Following the proof of Theorem 1.2, we need to check that $ is a contraction on 
X^^^([0,T]), the closed convex set of C*^([0, T]; Co(R''~)) consisting of functions w{t,r) such 
that w{t,r) + e~^vo{r) G C^{[0,T], X^). First we estimate \^{t,r,w{r))\, and write 

(1.29) mt,r,w)\ <^j^ dt^e'^-WAnw{h,r)\ + \v^{r)\ + e-'^\AnVQ{r)\) 

By Proposition a.2, A^ is a bounded operator on X^, and since vi G X^ by assumption, 
we have w G X° ; by the same remark as in the proof of Theorem 1.2, we conclude that 
r, w) G Then (1.23) shows that $ is a contraction on X^^„^([0, T]), when T > is 

small enough, and the proof goes as in Theorem 1.2. ^ 

Again, by the group property, we find V G C°(R+,X^). From Theorem 1.5 we can infer 
existence of limit points of the orbits : 

Corollary 1.6: With the same hypotheses as in Theorem 1.5, from any sequence tn — > +oo, 
we can extract a subsequence tn^ such that v{tnj,r) — > G C°(R"'') as j — > oo for the 
convergence on compact sets. 

Proof. The family u{t, r) = mp + v{t, r) is clearly bounded by 1, since m(t, x) = e*"^w(t, r), 
X = re*^ solves the evolution equation corresponding to (1.10). Consider next w{t,r) as 
in the proof of Proposition 1.4, (1.23) shows that ^(A„(w(s,r) +fi(r) + e~'^^„vo(r)) is 
bounded uniformly in (t, r), so by integration of (1.25), drw{t,r) is uniformly bounded, and 
so \drv{t, r)\ < C. It follows that the family v{t, r) is also equicontinuous, and the conclusion 
follows from Ascoli-Arzela theorem. X 

Let us extend once more our previous considerations. We enrich the structure of our 
Banach space X^ by requiring some asymptotic behavior near oo. Consider indeed the set 
of functions u G Xf, such that ru{r) has a limit as r — >^ oo, and if (,{u) = WvHr^oo ru{r), 
the function r{ru{r) — £{u)) is bounded. It is easy to see that F° is a separable Banach space, 
with norm 

(1.30) \\u; Y^W = sup^£[o,i[ Hr)\ + sup^g[i^+<^[ \ru{r)\ + sup^g[i^+<^[ \r{ru{r) - £(«)) | 
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Similarly, consider the set Y2 of functions v e , such that rv{r) has a limit £o('J^) as r — > 00, 
r^v(r) — rlQ{v) has a limit £i{v) as r — > 00, and the function r^/'^{r'^v{r) — r£o(v) — ii(v)) is 
bounded. It is easy to see that is also a separable Banach space, with norm 
(1.31) 

\\v;Y^\\ = sup^g[o,i[l'f^('')l + sup^e[i,+oo[k^^(OI + sup^e[i,+oo[k^^^(''^«(0 -''^o(^^) -4(^^))| 

By 1^°, A; = 1,2 we denote also as before the closed convex subspace of Y^ consisting of 
functions v such that v{r) + mp>0, and yJ([0, T]) = C°{[0, T],Y^). 

Proposition 1.7: With the assumptions of Theorem 1.5 then for all T > small enough, 
and all vq G Y^, k — 1,2, there is a unique v e ^^''([O, T]) such that u — v + rrifs verifies (1.11) 
with u|t=o — Vq + mp. 

The proof goes along the same steps as this of Theorem 1.5, but this time we need also 
Proposition a.3 to control the last term in (1.30) or (1.31) for ||t)(t, •); F^H. Note that we do 
not expect asymptotics beyond this order, on account of our accuracy in estimating An- 

Whatever the class or Y^ to which the initial datum does belong, we cannot ensure 
that V* itself belongs to this set, nor even to Co(R"''). Fortunately the main properties we 
shall use don't hinge upon v* itself (except for Theorem 0.2, where we are lead to make an 
hypothesis about the short range of v*, ) but are rather inherited from those of the finite 
time evolution v{t,r), t > 0. For this reason we shall also call a limiting orbit of the flow a 
X- ghost, stressing that it proceeds from an initial datum in X. Actually the phase picture 
for the whole dynamics may look quite complicated, unless we could prove uniqueness of the 
limiting orbits of (1.11), at least for a given uq. 

Now we extend the gradient-flow dynamics to those positive functions, which are merely 
continuous on R+ and bounded by 1. Indeed we shall eventually take a X^-ghost u* as a 
new initial datum, but use u*{t, r) only with r in a compact set. We have : 

Proposition 1.8: Let as above be positivity preserving, but not necessarily with the 
asymptotic property. Then for all T > small enough, and all uq G C'^(R"'"), ||wo||oo < 1, 
there is a unique u G C*^(R+)([0, T]) such that u{t, r) verifles (1.11) with ^1^=0 = ^o- 

The proof is omitted, for it goes as in Theorem 1.5, without the additional requirements 
on asymptotics at inflnity. Iterates of $(t,r, w) to the flxed point converge uniformly for r 
on every compact of R"*". Again, by the group property, we flnd u G C'^(R"'"; C*^(R"'")). 

Next we consider various comparison theorems. As in the scalar case [Pr] , the comparison 
theorem or maximum principle shows very useful in our situation. To flx the ideas, we state 
it for the full dynamics. Recall that = + mp,v^ G C^(iL'^,X^) is a supersolution 
of the Cauchy problem (1.11) with initial datum Uq G X^ if ||«"''(t, •)||cxo < 1, '1*0 > uq and 
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verifies 

^ > + f{PAr,U+) 

We define analogously a subsolution u~ . Because An is positivity improving, for all t > 0, 
u~{t^r) < u{t,r) < u'^{t,r). The first consequence of the maximum principle is that u 
doesn't increase beyond if this holds for uq. More precisely, we can show the following : 

Proposition 1.9: Let u{t,r) = v{t,r) + 771(3, v E C°(R"'",X^) be the solution of (1.11) with 
initial datum «o(r) satisfying |uo(^)| < A* < 1 for some > m^. Then \u{t,r)\ < jji and for 
aU t > 0. 

The easiest way of proving this is to consider m(t, r) = e^^^uit, r), and the argument goes 
as in as in [El-BRo, Proposition 3.3]. This holds equally for the partial dynamics. Another 
consequence of the maximum principle is monotonicity : 

Proposition 1.10: Let An be positivity preserving. If the initial datum tto(r) in (1.11) is 
an increasing function of r, then the same holds of tt(t, r) for all t > 0. Assume moreover An 
be positivity improving. If tto(r) is strictly increasing, then the same holds of u{t,r) for all 
t > 0. 

Proof. Let a > 0, and consider the new dynamics on r e]a, +oo[ given by = —Ua + 
f{PAnUa), Ua{0,r) = Ua, where Ua{r) = uo{r — a) is the translate of uq- Since Ua{r) < uo{r), 
we see that tt is a subsolution of ^ = —u + f{j3Anu), so u{t, r — a) < u{t, r) for all t > 0. X. 

Note this implies the former result if we take A = mp in Proposition 1.9. Of course these 
properties extend, by continuity, to the limit points u*{r). 

When considering the spectral problem in Sect. 3, we shall need to know that the limiting 
orbits belong also to C^{R.'^). Taking the r-derivative of u in (1.11) involves the derivative 
of AnU which we can compute using (1.24), but it is very hard to give estimates on Bn along 
the lines of Appendix A, since Bn doesn't enjoy the asymptotic property and so forth. Since 
we will eventually take as an interaction the function J of Example 1, we restrict to this case, 
which is much simpler because the closed form of Weber second exponential integral involves 
just one (modified) Bessel function /„, instead of the correlations Jn{r p) Jn{r' p) . 

As in (1.18), we consider the space Xl of functions v e C^{Ti'^) fl such that 
supj,^[i \r'''^^v'{r)\ < oo. The weight r'^'^^ is chosen in such a way that we can take the 
derivative of the asymptotics of v, and is consistent with the complete asymptotics (1.13). It 
is easy to see that Xl is a separable Banach space, with norm 

(1.34) \\v;Xl\\ = ||z;;XO||+sup,g[o,i]l^'WI + suPr6[i,+oo[k''+'^'(OI 

For /c = 2, we have X^, C (the usual Sobolev space. ) Let also as before be the closed 
convex subset of Xl consisting of functions v e Co(R''"), such that v{r) +m/3 > 0, and define 
Xl{[0,T]) = C%[0,T],Xl). We have: 
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Theorem 1.11: With the assumptions above, for all T > small enough, and all vq G X^, 
k = 1,2, there is a unique v e X^([0, T]) such that u = v + mj3 verifies (1.11) with u\t=o = 
Vq + mp. Moreover, from any sequence tn — > +oo, we can extract a subsequence t^j such 
that v{tny.,r) — > V* e C^{R^') as j — > oo for the convergence on compact sets. 

Sketch of the proof. We argue as in Theorem 1.5, showing that for small t > 0, $(t,r, ■) is a 
contraction on X], . This follows from the fact that (in the particular case where J is a 
Gaussian) is a bounded operator on X\, as can be shown by using the asymptotics at infinity 
of In as in Proposition a.2 : for large rr' we replace using (1.13), exp[— (r^+r'^)/47r]/n(rr'/27r) 
by {rr')~^/'^e~^'^~'^ g^j^j ^.g^y standard gaussian integral arguments. 

The last part of the Theorem follows as in Corollary 1.6 from the equicontinuity of 
the second derivatives of Anu{t,r), as we can check by iterating (1.22), and Ascoli-Arzela 
theorem. ^. 

Note that, since we proceed by extraction of subsequences, the -ghost we obtain 
from the sequence v{tn., r) as a limiting orbit in the C^-topology may not coincide with the 
corresponding X^-ghost in the C'^-topology obtained in Corollary 1.6, even with the same 
initial datum. But the -ghosts also enjoy the monotony property as in Proposition 1.10. 

b) The partial dynamics, and the Barrier Lemma. 

In order to cope with divergent integrals we need also study the dynamics in some finite 
boxes A = {|a;| < A}, for which A — > oo. Outside A, the magnetization m is frozen to a 
configuration m^^c which acts as a boundary condition for the evolution inside A. We follow 
again closely the main steps of [Pr] . Define the free energy with boundary condition as 

(1.40) .^("7,a|"7-a=) = .^a("^a) H — dx dx' J{x — x')\mA{x) — mA'={x')\'^ 

2 J A Ja'= 

where 



(1.41) .^a(?t^a) = 7 I dx I dx' J{x — x')\'mA{x) — mA{x')\^ + I dxfpimA 

4 Ja Ja Ja 



(^)) 



[Contrary to Appendix C, we have removed from the second integral the term f{l3mp), 
which amounts to shift J-'a{'>^a) from a constant term, so long A is kept fixed. ] Since 
J J = 1, J-'{mA\mA'^) < oo for all A. We take variations of JF inside the sector (e„), i.e. 
among all radially symmetric functions of the form m{x) = u{r)e^'^^ . To simplify notations, 
remove subscript n from and define operators Ax and by their kernel A\{r,r') = 
x{r' < A)A(r, r') and A^(r, r') = x(r' > X)A{r,r') respectively (here x denotes the (sharp) 
charateristic function. ) We denote also hj u = ux Q) uj^ function equal to ux on [0, A] 
and outside, and write Au = {Ax ® A'^){ux © u'^) = Axux + A'^u'^. We shall use twidled 
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u and V to refer to partial dynamics. It follows from (1.40) that Euler-Lagrange equation for 
J^{-\mAc) restricted to the sector (e^) is given by 

(1.42) F^{u) = -Au + ^I'{ux) = 
which is equivalent to 

(1.43) F^{u) = -ux + f{(3Au) = 

So we define the partial dynamics by T^uq = u\{t, ■) ® uj^ where ux{t, •) solves 

(1.44) ^ = -ux + f{pAxux + PAlul), ux{0,r) = uo{r), 0<r<A 

with the same initial condition uo{r) as for the full dynamics (1.11). The integrated form of 

(1.44) is again 

(1.45) ux{t,r)^e-'uo{r)+ [ dse'-' f{(3AxUx{s,r) + (3A'xul{s,r)), 0<r<A 

Existence and uniqueness for (1.45), expressed in terms ofv = vx®v1 for the decomposition 
u = {mp + vx) © [171/3 + vl) follows as before, and vx G C°(R+; W^{A) n C^(A)) [so long as 
we are concerned in the first r-derivative. ] Compactness of the orbits results also from the 
uniform boundedness of vxit^r) and drVx{t,r). Namely, given any sequence tn oo, there 
IS uX e Ci(A) such that limfc^oo 

u = u* for a subsequence — > oo. 
Note that Propositions 1.9 and 1.10 also apply to partial dynamics (1.44). Decay prop- 
erties of as A — > oo will follow from the "Barrier Lemma" (in the terminology of [Pr], ) 
which is an essential tool in our analysis ; this compares the full dynamics as in (1.11) 
with the partial dynamics T(^. 

Theorem 1.12: Let u{t, r) and u{t,r) solve respectively the full and partial dynamics (1.11) 
and (1.44), with same initial condition uq G X2 . Then for all T > 0, u{t, r) — u{t, r) and 
dru{t, r) — druit, r) — > uniformly for < r < A and < t < T, as A — > 00. [the r-derivative 
of u{t, r) being understood almost everywhere, namely outside r = A.] 

Proof: We proceed somewhat as in the proof of [El-BoRo,Proposition 3.5]. Recall that the 
restriction of to [0, A] is given by 

(1.50) ux{t,r) = e-'uo{r)+ [ dhe-^'-*^^ f{pAxUx + pAlul){tur), r<X 

Jo 

where u'^ is independent of t and takes for instance the initial value Uq. [Since we are not 
interested in uniqueness properties of the limiting orbits, we could assume = m^j. ] Denote 
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by u(t, r) the solution of (1.11) on the full space and by tt(t, •) = T^uq = ux(t, •) ® tt^. So u 
is piecewise continuous. We have 

(1.51) {u-^{t,r) = x{r<X) f dt^e-^'-''^ gxM{ti,r) + x{r > X){u - ul){t,r) 

Jo 

with gx{u, u) = fipAu) - fipAxux + /3^>^). Since < /' < 1/2, we get 

(1.52) \9x{u,u){t,r)\ < ^\A{u - ul){t,r)\ = ^\Ax{u - ux){t,r)+Al{u - ul){t,r)\ 

all r < A. Proposition 1.3 shows that if Uq — mp G X^, then u — rufj G X]^([0,T]) for all 
T > 0, and the RHS of (1.52) is weU defined. Applying A = Ax®A\io (1.51) we get in turn 

A{u-u){ti,r) = Ax / dt2e'^-''gx{u,u){t2,r)^Al{u-u'x){ti,r) 

Jo 

Using (1.52) and the fact that Ax is positivity preserving, gives the estimate 

\A{u-u)\{ti,r) < J^' dt2e''-''Ax\A{u-u)\{t2,r) + \Al{u-ux)\{ti,r) 

Inserting into (1.51) we get 

\u-u\{t,r) <x{r < X){-)^ [ dtie'^'-'^^ [ ' dt2e-^''-''^ Ax\A{u - u)\{t2,r) 

2 Jo Jo 

+ x{r<X)^jyhe-('-''^\Al{u-ul)\{h,r) + x{r>\)\u-ux\{t,r) 

Let T^'^^u{t) = Jq dti Jq^ dt2 ■ ■ ■ Jq''~^ dtku{tk) denote the /c-fold integral of u, this formula 
can be carried over by induction as 

\u - u\{t, r) < x(r < X)e-' E(f )'^^''^ (e(-)yir Va(« - «a)I) (^, 
(1.54) 3=1 ^ 

+ {^^fT^k) (e(-)^J-i|^(^ _ u)\){t, t)] + x(r >X)\u- ul\{t, r) 

We first need an estimate on A^^{u — u'j^){t, r) for r < A. We proceed as in Lemma c.l, using 
(c.6) 



(1.55) 



Al{u-ul){t,r)^{n + l) / dr'{u - ul){t,r') / dpTJ{p)Jn{rp)Jn+i{r'p) 

J\ Jo 

+ / dr'{u-ul){t,r') / dppJ^J{p)Jr,{rp)—Jn+i{r'p) 
Jx Jo "P 
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We split dpJ^J{p)Jn{rp)Jn+i{r'p) into 2 parts, integrating respectively on [0, 1/A] and 
[1/A, oo[. For the first one, we make use of the bounds \ Jn{rp)\ < C{rp)'^, Jn+i{r'p) < C, 
for the second part, of the bounds |J„+i(r'p)| < C{r' p)~^^'^, Jn{fp) < C, and of the rapid 
decrease of TJ{p). Altogether, we get with a new constant C > : 

(1.56) I dprj{p)UTp)Jn+i{r'p)\ < c{^{{^T ^ ^) + 

Next we integrate by parts Jg°° dpp!FJ{p)Jn{rp)^Jn+i{f'p)i the 2 first terms can be bounded 
as before, for the third term j^dppr!FJ{p)J'^{rp)Jn+i{r'p) we split again according to 
[0, 1/A], [1/A, oo[, for the first integral we use rpJ'^{rp)\ < C{{rp)'^ V ^/fp)■, which gives the 
bound C';x((l)"^ ^ V^^' ^ ^^^^ I'^nl'^P)! ^ ^i which gives the bound 

Cr/Vr'. Thus 

(1.57) I / dppJ^Jip)Urp)-^JnMr'p)\<C{{{{{)''yJ-) + -^) 
Jo clp A A \ p Vr' 

On the other hand, since the initial datum uq belongs to X2, given T > 0, by Theorem 1.5 
there is Ct > (also depending continuously on T) such that for alH < T, \u — tt^|(t, r) < 
Cr(l + r)~^. Inserting this estimate into (1.55), using (1.56) and (1.57) we get by integration, 
for a new constant Ct > 

(1.58) \Al(u-u'i)it,r)\<CT-^, r<X 

Next we notice that since e*"'^^„tt(r) = J*(e*"''tt), and J J = 1, we get j^A'^^K'^) < sup[Q ;)^] \u\ 
for all r > 0. It follows then from (1.58) that 

A^-'lAliu - ul)\{t, r) < Ct/^, J = 1, 2, • • • 

Now we have again \A{u — u){t, r)| < 1 for all t,r >, and A'^~^\A{u — u){t, r)| < 1 for all k. 
Performing the successive integrations we find that the series in (1.54) is uniformly convergent 
for t in compact sets, so we can write, for < r < A, 

\u - u\{t, r) < e-' f:(f )^T(i) (e^-)^- ^ |A^(« - <)|) {t, r) < Cre^W'-'^ 1= 

j = l V 

which proves the first estimate on u{t, r) — u{t, r). 

We cannot repeat this argument for the r-derivative, since dr doesn't commute with Ax 
[while partial derivatives dx^ commute with convolution, ] so we proceed indirectly. Take again 
r-derivative of (1.25), and iterate (1.22) once more to get \d^Anu{t,r)\ < || J||i ||u||oo- By 
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integration in the s variable, we find that d^w{t, r), and consequently d^u{t, r) are uniformly 
bounded for t in compact sets. The same conclusion holds for the partial dynamics, showing 
that \d^(u{t, r) —u{t, r)) | < Const, uniformly for t in compact sets. Now we use the following 
well-known interpolation inequality : for any bounded intervals Ki (ZC. K2, there is C > 
such that for all f E we have 

sup^^ l/f < Csup^^ |/|(sup^^ I/I +sup^, |/"l) 

Applying this to /(r) = u{t,r) — u{t,r) easily yields the required uniform estimate on 
dru{t, r) — dru{t, r). The theorem is proved. X 

2. Existence of a radially symmetric minimizer. 

In this Section we study some continuity properties of the free energy in finite or infinite 
volume ; continuity properties of the renormalized free energy are investigated in Appendix 
B. Next, following [Pr], we identify the limit points u*{r) in the box A by using the excess 
free energy functional .F(mA|mAc), then we identify the limit points u*{r) on the half-line, 
by using the renormalized energy J-ren{fn) in the full space, and eventually prove Theorem 
0.1. 

a) Continuity properties for the free energy. 

We start with general remarks on continuity properties for the free energy, which are 

quite close to the case of configurations valued in [—1,1] as in [Pr], [AlBe], The only 

difference is that we need to assume |m| 7^ 1, since the entropy function I{m) is unbounded 
as |m| — > 1. Consider the functional 

(2.1) T{m) = - I dx I dyJ{x-y)\m{x)-m{y)\^+ ! dxfp{\m{x)\) 

defined on the set E = L'^{B?] ^2(0, 1)) and valued in R+ U {+00}. 

Proposition 2.1: T is lower semicontinuous on E with respect to convergence almost ev- 
erywhere. 

Proof. Let m E E, and & E he a sequence converging a.e. to m. Define 

9n{x) = ^ / dyJ{x -y)\mnix) -mn{y)f + ff3i\mn{x)\) 

and g{x) similarly, with mn{x) replaced by m{x). Expanding the square and using the 
normalization of J in we get 
I 

9n{x) - g{x) = -(|m(a;)p - |m„(a;)p) + f^{\mn{x)\) - fp{\m{x)\) 

— - Ke{mn{x) — m{x))J * m{x) + - Rem„(a;) J * (m^ — m){x) + -:J * (l^nl^ ~ l™l^) (^) 
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The first term tends to a.e., so does /g(|m„(a;)|) — fp{\m{x)\) since is continuous. Since 
m e L°° and J & L^, J * m{x) is uniformly continuous and bounded, so again {mn{x) — 
m{x))J * rri{x) tends to a.e.. In the same way, by the dominated convergence theorem, 
J * (m^ — m)(a;), tends to locally uniformly in x, and this argument equally applies to 
J * (|mnp — |mp)(a;), proving that the last 2 terms in the decomposition of gn{x) — g{x) 
above tend to a.e.. Then Fatou lemma shows that 



(2.2) liminfn^oo .^(^T^n) > J^irn) 

which proves the Proposition. X 

Thus proving the existence of a minimizer for T amounts to extract from every mini- 
mizing sequence a subsequence converging a.e.. The next result concerns weak lower semi- 
continuity of the free energy in finite volume A as is defined in (1.40) and (1-41). To simplify 
the notations, we set m = rriA, and m*^ = rriAc. Here m'^ = rriAc E E is fixed. 

Proposition 2.2: If G E converges weakly to m in Lp{A), 1 < p < oo, then 

(2.3) liminf^^oo •^("^nl"^'^) ^ T{m\m^) 

while, if |mn(a;)| < ji < 1 and nin — > m a.e. in A, then J^{mn\m'') — > J^(m|m'^). 
Proof. Following [Pr] we write J^{m\m^) = F{m\m^) + R{rrf) where 

F{m\rrf) — — / dxl{m{x)) Re f dx dyJ{x — y)m{x)m{y) 

P Ja 2 J a 

(2.4) — Re / dx dyJ{x — y)m{x)m'^{y) 

J A J A'' 

R{m^) — - I dx I dyJ{x — y)\rrf{y)\^ 

So R{m^) is just a constant. Because G E converges weakly to m in L^(A), j^dyJ{x — 
y)mn{y) tends to j^dyJ{x — y)m{y) as n ^ oo for any a; G A. Moreover Ii{m) = 
I Re dx dy J{x — y)m{x)m{y) is weakly continuous on Lp{K) ; namely, if m„ converges 
weakly to m, then mn{x)mn{y) converge to m{x)m{y) weakly* in L°°(A x A), and since 
J{x — y) belongs to L^(A x A), then Xi(m^) Xi(m). We get the same conclusion for the 
term X2(m) = dx f^^ dyJ{x — y)mn{x)m'^{y), so that (2.3) follows from the convexity of 
the entropy function I{m) = I{\m\). 

If nin rn a.e. in A, then by the dominated convergence theorem, the last 2 terms in 
F{mn\m^) tend to the corresponding ones with m instead of nin- The same conclusion holds 
for the first term provided \mn{x)\ < // < 1 in A. X 
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Following [AlBe,Remark 4.8] we can use the first part of Proposition 2.2 to prove existence 
of a minimizer for J^{-\m^) subject to some weakly closed constraint, since J^{-\m'^) itself is 
not coercive. We can also extend Proposition 2.2 to infinite volume for as in (2.1), to show 
that E is weak* compact in L°°(R)^, and is weak* lower semi-continuous on E. Thus we 
can prove again existence of a minimizer for T subject to some weakly closed constraint ; but 
it turns out that a topological constraint such as the degree of m at infinity is not weakly 
closed, and thus we shall proceed another way. 

b) Free energy dissipation for the partial dynamics. 

This paragraph is a first step towards Theorem 0.1. The key property of the solution of 
(1.11) or (1.44) is that its energy decreases with time. To start with, we consider the case of 
finite volume. Given A we write again, emphasizing the dependence on 9 in the sector (e„), the 
solution of the partial dynamics T-^ with initial value e*"'^'Uo(r) as m{t,x) = T/^{e^^^uo{r)), 
using the notation m ~ m\ © = e^'^^u = e^'^^{u\ © -u^). Sometimes we omit to write 
gin6i_ -yy^g define the free energy dissipation rate of u{t^ r) as 

(2.8) X^(2)(t, r)^^£ dr r{-/3Au + T'{u^)) {ux - f{/3Au)) 

It is easy to see that J^{mA\m\c) is a Lyapunov function for Eqn. (1.44), i.e. X^{u) > and 

(2.9) {(T^fh) A\mAo) - T{mA\m Ac) = - [ dsl^{u){s,r) 

Jo 

with I^{u) = iff w verifies (1.44). If the initial datum u is bounded below from 1, so is ux 
because of Proposition 1.10, and I'{u\) < oo everywhere. This shows that X^(tt)(t, r) < oo. 
We have the following 

Theorem 2.3: Let the initial datum m be such that J^{mA\mAc) < +oo. Then every limit 
point u*{r) of T^u satisfies Euler-Lagrange Eqn. (1.42)-(1.43) and 

(2.10) J^(e^"%^(r)|mAc) < J^(mA|mAc) 

Proof. We follow an argument of [Pr, Sect. 4. 2. 5], essentially due to [FiMc-L]. Since u* is a 
limit point of T^uq, there is a sequence tn oo such that lim^^oo \\T-^uo — tl*||oo = 0. If 
u* does not satisfy (1.42), then X^(u*) > 0. We start from u* as a new initial datum for 
the evolution T^. Because t i— T^{u*) is continuous as a map C'^(A), and is 

continuous (hence l.s.c) on C°(A), for t > sufficiently small we have : I^{T^u*) > and 
hence the free energy dissipation in [0, 1] verifies 

(2.11) D^(u*) = [ dtX^{T^u*) > 

Jo 
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We shall show in a while that this implies 

(2.12) limt^oo / dsI^{TSo) = +00 

which contradicts hypothesis J-'{m\\m\c) < 00 by inequality < in (2.9) : 

f dsX^{T^uo) < r{mA\mKo) - .F((T/-e*"%o)A|^^Ac) < T{mA\mAo) < 00 
Jo 

So for r < A we have, —j3Au* + /'(w^) = 0, or equivalently u\ — f{(3Au*) = 0, and the 
limiting orbit u* verifies Euler-Lagrange equation (1.42). Then (2.10) easily follows from the 
uniform convergence on compact sets, as t — > 00, of {T^e^'^^uo)x towards e*"^tt^. 

Now we need to show that (2.11) implies (2.12). Because of the continuous dependence 
of the initial data, we have again lim^^oo sup^<^ \\T^u — T^{T^uo)\\oo — 0. Since we work in 
the finite volume A, and T^uq is bounded below from 1, this implies lim^^oo D^{T^uo) = 
D^{u*) > 0. So there is iV G N such that D^{T^^uo) > 5 > for n > iV. Without loss 
of generality (since — > 00) we can assume, after possibly extracting a subsequence, that 
N = 1, to = 1 and tj — > 1 for all j > 1, so ti > 2 and by the group property we have 

dtX^{T^uo) = / dtX^{T^T^uo) > D^{T^uo) > S 
Jo 

By induction we get when t — > 00 : 

/ dsX^{TSo) > ( / ' + / ' + • • ■)dsX^{T^^uo) >S + S + ^00 

Jo Jl Jti 

which proves (2.12). X- 

c) Free energy dissipation in infinite volume, and proof of Theorem 0.1. 

We extend here the results of Paragraph b) to full dynamics. Let m E E he the initial 
condition for the full dynamics, such that J-renim) < 00. More precisely, we choose m{x) = 
e*'^^wo(r), uo = vo + mis e W, vq G with notations of Theorems b.6 and 1.5. For A > 0, 
denote by m = rriA © wa^, and u = ux + uj^ the corresponding radial parts. Here we shall let 
A — > 00, for fixed t. Using that {T^m)\^c = mAc, Theorem b.6 and Remark b.7 easily show 
that we can rewrite (2.9) as 



{u){s,r) 



(2.15) J^ren(7fm))--^ren(m)<- / dsX 

Jo 

By the Barrier Lemma, Theorem 1.12, we have \\T^rh — Ttm; L°°{A)\\ = supj.^[Q \u{t,r) — 
u{t, r)\ — > 0, and \\drT^m — drTtm; L°°(A)|| = sup^^^.A] \dru{t, r) — dru{t, r)\ ^ 0, as A — > 00 
uniformly for t E [0,T]. By [lower-semi] continuity of J^rem Theorem b.8 

(2.16) liminfA^oo [^ren(Tfm) - J^ren{m)] > Tren{Ttm) - J^ren{m) 
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On the other hand, Fatou Lemma shows that 

hminfA^oo / dsX^{u) > I dsI^{TsUo) 
Jo Jo 

where is defined as in (2.8) by integrating over r G [0, oo[ and takes its values in [0, +oo]. 
We rewrite this inequahty as : 

hmsup;(^^QQ— / dsl^(u) < — / dsI^{TsUo) 
Jo Jo 

so by (2.15) 

(2.17) Tren{Ttm) - Trenim) < - f dsI^{TsUo) 

Jo 

Now let m*{x) = e'"^tt*(r) be a limit point of T^m, u* G C°(R"'") as in Corollary 1.6, namely 
suppose there is a sequence —>■ oo such that for any A > 0, lim^^oo sup^^;,^ \Ti^uo{r) — 
u*{r)\ = 0. So this time we let t ^ oo, for fixed A. If u* does not satisfy (1.10), then there 
is A* > so that I^* (u*) > 0. 

We start from u* as a new initial datum for the evolution Tt. Note that u* has a 
priori no decay at infinity, but we can consider Ttu* on r G [0, A*] by Proposition 1.8, with 
convergence on compact sets. Because t i— > Tt{u*) is continuous as a map R"*" — * C^{A*), and 

is continuous (hence l.s.c) on C'^(A*), for t > sufficiently small we have : (Ttu*) > 
and hence the free energy dissipation for (t, r) G [0, 1] x [0, A*] verifies 

(2.19) D*{u*)= [ dtl^* {Ttu*) > 

Jo 

As in the proof of Proposition 2.3, we shall show that this implies 

(2.20) limt^oo / dsl^* (TsUq) = +oo 

Jo 

On the other hand, inequality (2.17) yields dsl^* (TsUq) < Jq dsI(TsUo) < J^reni'm) — 
^Ten{Ttm). So by the discussion at the end of Appendix B, there is C > such that 
•^ren(Ttm) > — C for alH > large enough, and we see that (2.20) contradicts the hypothesis 

y^venim) < oo. 

Now we prove that (2.19) implies (2.20). Since Tt^uo{r) — > u* uniformly for (t, r) G 
[0, 1] X [0, A*] as n ^ oo, and the flow Tf depends continuously on the initial data, we see 
that Tt{Tt^uo){r) - Ttu*{r) ^ as n ^ oo, uniformly for {t, r) G [0, 1] x [0, A*], and by Fatou 
Lemma, 

lim„_oo^*(Tt„wo) > / dtX^* {TtU*) > {) 

Jo 
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So we may assume that D*{Tt^uo) > 5 > for all n large enough, and the proof goes exactly 
as in Proposition 2.3. This shows (2.20) and brings the proof of Theorem 0.1 to an end. X 

3. Linear stability. 

To start with, we recall from [OvSil] some well known facts concerning symmetry 
breaking of Eqn. (1.8) or (1.9). The symmetry group G for these equations is given by 
G = X 0+(2) X U{1) X r, where acts as translations m{x) m{x -h), h e B?, 0'^{2) 
as rotations m{x) — > m{R~^x)), R e 0'^{2), U{1) as gauge transformation m{x) — > Xm{x), 
A e t/(l) S^, and F as "conjugation of charge" m{x) m{x). 

By the symmetry group Gm of a solution m, we mean the largest subgroup of G which 
leaves m fixed. Then the part of G broken by m is the coset G/Gm- If is a one-parameter 
sub-group of G, we say that H is preserved (resp. broken) by m if h{m) = m for all /i e 
(resp. h{m) ^ m for all h E H,h^ Id. ) 

The subgroup of translations is never preserved by m, unless m is a constant. The sym- 
metry group of a radially symmetric solution, i.e. m{x) = e*"'^tt(r) is the discrete subgroup 
of 0^^^(2) C 0'^{2) of rotations by the angles 2k7r/n,k e Z. Thus, m breaks the trans- 
lation group, the rotation subgroup 0'^{2)/0^^^{2), and the charge group. The symmetry 
group for an equation of the form F{m) = allows to find elements in the kernel of its 
linearization around some point m, i.e. solutions of {dF{m),$,) = 0. Namely, let m be a 
solution of F{m) = 0, breaking a one parameter subgroup g{s) e G (the symmetry group of 
this equation). Let r be the generator of g{s). Then ^ = rm solves the linearized equation 



a) Linearisation of Tren around a radially symmetric solution. 

Next we examine the Hessian of J^^en and look for relations between the gradients of 
Fi and Fq defined in (1.8), (1.9). We take advantage of the fact that Fj, j = 0, 1, are real 
analytic functions, real for real m, to write Fj{m) instead of Fj{m,fn), and we will denote 
also Fjim) = Fj{m,rn) = Fj{m,m), ^^{tti) = ^^(m), etc . . .Following [OvSi], we express 
the Hessian of Tren in these complex variables as 



We have Hess.7ren("7.) = VmFi{m) (the gradient in the real sense) and similarly, we define 
VFo(m). Introduce the notations z = (3J * m, (j)Q{z) = /(|2;|) = , 0i(m) = /'(|m|)|^, so that 
00 and 01 are inverse from each other. 

Lemma 3.1: With the notations above 



(dF(m),e) = 0. 



(3.1) 




(3.2) 



Id +/3V^0o(-2) J * • = V™Fo(m) 
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and in particular (^PV z(l>o{z)) ^ = VmFi{m) + J*. 

Proof. Again, we write 0i(m) for 0i(m, m), (j)i{fn) for (^i{m,m), and similarly for (f)Q. The 
"upper- left" matrix element of V0o(-2) VFi(m) is given by 

On the other hand, differentiating the identity 0o o 0i = Id we get 

1 a iM M \( \ d(f)o d(j)i d(po .d(j)i 

1 = 9^(00 o 0x)(m) = + ^(-) 5=M 

which leads to Pai = 1 — f3^^{z)J * ■. All other matrix elements can be handled of this sort, 
using drn{(l>o ° 0i)("^) = 0, and two similar identities, obtained after permuting z with 'z, m 
with m. So we proved 

(3.3) -/3V0o(^)VFi(m) = -Id+/3 ^ _ g _ J* • 

and the Lemma easily follows. X 

A direct computation also shows ^^{z) = K^^j^ + /'(l-^l))- Using recursion formulas 
between the derivatives of the modified Bessel functions Iq and /i, we find for / = Ii/Iq : 

f(t) = i-p{t)-i^,so 

(3.4) ^-^{z) = \{l-f\\z\)) 
Similarly 

(3.6) ^(.) = ±(l_A|,|)_2/(M) 



SO that setting 6(|^|) = 1 - f{\z\) - 24^, we have 



(3.7) VAo{z) = - 



1 (i-f\\z\) mz\) 



z 

2/ 



^^h{\z\) 1-P{\z\) 
It follows that V z(l>o{^) is hermitean, and 

(3.8) det V0o(^) = a{\z\) = /'(|^|)4t^ 

We notice that 1 - P{\z\) > 0, and 6(1^1) = f'{\z\) - < with equality only at ;2 = 0, 

since /(O) = and / is strictly concave on [0, +oo[. Now by Lemma 3.1 



(3.9) VFo(m) = -Id + ^ ( ^ , ^"l^Zn ) J* 



P fl~fH\z\) §b{\z\) 
lb{\z\} 1-P{\z\) 
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Let now m be a radially symmetric solution of dUFy-en = 0, as in Theorem 0. If m(a;) = 
e'^'^uir) solves (1.8), then f(\z\) = uir), and we have 6(1^1) = hir) = 1 - u^ir) - 2..^^, 

I'ou{r) 

a{r) = a{\z\) = f o f^^ o u{r) j-^^^^) (here we have used the relation \m\ = f{\z\) to denote, 
somewhat incorrectly, a{\z\) by a(r), h{\z\) by h{r). ) Then Lemma 3.1 shows again 

(3.10) HessJFren(m) = V^Fi(m) = -^B{r,e) + J* • 



with 



(3.11) B{r,e) 



1 f l-u^{r) -e2^^^6(r 



o(r) V-e~^*"^&(0 l-u^r) 
(cf. [OvSil, formula (7.7)] for Ginzburg-Landau equation. ) As in [OvSi] we have the 

Lemma 3.2: VFi is symmetric for the scalar product r]) = Re J rj^. In other words. 
Re / f]WFi{^) = Re JVi{r])^, where VFi(0 is a shorthand for {VFi{x), (C,C))- 

Now we make a Fourier analysis of V-Fi. Taking polar coordinates as above, write 

^ = ^k{r)e^'^^ , ^fe G C, and expand J * ■ in terms of A^. To account for the phase factors 

fcez 

g±2in0 (3.11) we make a shift of indices and introduce the mapping 

^ ■■ (C,0 = ((6,4))fcgz = ((^fe'?2n-fe))fcez 

which is unitary if the target space is endowed with the inner product 



(3.13) {lv)-MLvn)+ReJ2{L^' Yf-"^' ) 

l^y^ \S2n-fc/ \V2n-k/ 



) 



We then define the real-linear operator VmFi on vector-valued functions ^ by Vm-^i = 
TrVmFiTT*. As in [OvSi] we can easily show that operator VmFi is block-diagonal of the form 



(V„F„f)=((V„Ff',( 

\?2n-A;/ 



where 

(3.14) V^Ff^=| ^^'^^g^ ^ ^7^^:^ 0. |=^A 



^fe 2^3^(1 w^(r)) 

2^a(l) ^2n-fc - 2^9^(1 - W^(?')) 



So we have reduced the problem of linear stability around the radially symmetric solution 
m{x) = e*"'^w(r) to the study of the spectrum of the family of operators {Lk)k>n- By the 
discussion at the beginning of this Section, we already know some zero modes, corresponding 
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to breaking of the rotation, gauge or translation symmetry. However, these functions are not 
in L2. 



b) Spectral properties of L^. 

Conjugating with ^ ~ i i ^ ' thai L„ is unitary equivalent to 




(3.15) 



Let V{r) = ^ > 0, W{r) = > 0. These are smooth functions on R+. Since 

AnU — ^I'{u) — U = (q) is a solution of LnU = (the zero mode is due to breaking the 
gauge group) but of course u ^ L'^. First we look at the spectrum of A^. 

Lemma 3.3: If J^J > and J > (i.e. J is non negative definite in the sense of [PrTo], ) 
then ]0, l[c ac{An) C a{An) C [0, 1]. Moreover Ker^„ = 0. 

Proof. RecaU An = HnJ^JHn, so (An - X)'^ = -jHn{l - ^y^Hn. So {An - A)-i is 
bounded in operator norm for A < 0, if TJ > 0. On the other hand, TJ{p) < 1 since 
J > and / J = 1, so when A > 1, (l - ^^Y^ is a Neuman series converging in 
operator norm. It follows that a{An) C [0,1]. Conversely, let A e]0, 1[. Since J(0) = 1 
and TJ{p) — > as p — > oo, there is pA > such that !FJ{p\) = A, so take v e L^(R+;r(ir) 
such that Hnv{p) = x{P\ ~ ^ < P < P\ + where £ > is small enough. It is clear that 
{An — X)~^v ^ L^, which shows that ^„ — A is not surjective. 

At last, let e L^(R+;rdr) such that Anip = 0. Then J^J{p)Hnip{p) = a.e. and this 
implies in turn that Hnip{p) = a.e. for we cannot have J- J{p) = on a set of positive 
measure since J is smooth and J > 0. So = and Ker A„ = 0. 4^ 

Consider now A^ — V{r), and A^ — W{r). We discuss first some properties of V and 
W . Convexity of the entropy function / on [0, 1] and /'(O) = imply that I"{u) > 0, and 

V'{u) = I"{u) ^ > 0. In the same way, it is easy to see that I"'{u) > 0. Recall also that 

/"(O) = 2, and mp satisfies the equation I'{mp) = /Srrip. Although this will not be needed, 
we mention for completeness that I'^mp) = {^^{/Smfs) — m^) 

On the other hand we learned from Sect. 2 that u increases from to m/3 on R"*". Alto- 
gether, this shows that V increases from 2//3 to 1 on r G [0, +oo[ for (3 > 2 (recall that there 
is nothing to prove when /3 < 2, ) and W increases from 2//3 to I"{mp). 

So in the sense of self-adjoint operators, we have An — 1 < ~ V{r) < An — 2/(3, and 
Lemma 3.3 easily shows that cr(A^) C [—1,1 — 2//3]. For convenience we shift the whole 
spectrum by 1 by changing V into V — V — 1, and consider now An — V. For a suitable 
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interaction J we shall prove that the spectrum of — V is purely continuous by using a 
positive commutator. This is actually a simple variant of Mourre estimates. More precisely, 
we assume that J(r) = (27r)^e~^^/*'^, or equivalently TJ{p) = e~'^^^ (see Example 1 of 
Sect.l. ) 

Lemma 3.4: With J as above we have [J * -jxdx + dxx] = —AttAJ, or in the (en)-sector, 
[A,, rdr + drr] = A'KH^{p^e--p")H^. 

Proof. We compute for a test function (/?, [J*-,xdx+dxx\tfi{x) = —2 f dxJ{x—y){x—y,'Vi^{y)). 
When J{x) = (27r)^e~^ Z^'^, —2J{x — y){x — y) = A-nJ'{x — y), so Green formula shows that 
[J * ■, xdx +dxx] — —4:TtAJ. Taking Fourier transform and restricting to the (en)-sector gives 
the Lemma. X 

Denote by D = rdr + drf the generator of dilations, and C„ = ATrHn{p'^e~'^f^)Hn. So 
Cn is a bounded, positive operator in the sense (and also positivity improving as we can 
show by differentiating (1.12) with respect to p at p = tt ; also Cn cuts off both low and high 
frequencies. These properties however, will not be used in the sequel .) 

On the other hand, [—V, D] = 2rV'{r), which is also a positive (and positivity improving) 
operator, so we get : 

(3.17) [An -V,D]=Cn + 2rV\r) 
From this and the energy identity 

(3.18) {{An - V, D]^\^) = {Cn^\^) + (2ry'(r)^|V) 

it follows that ij) e T^{D) = {■?/'£ L^(R+;rcir) : r'ip'{r) e L^} cannot be an eigenfunction of 
An — V, for otherwise it would be supported at r = (and the range of its frequencies would 
reduce to p = 0. ) But since D is not bounded on we shall first perform a regularization. 
This is simpler than usual (with the laplacian instead of a convolution operator) because An 
itself is bounded. Namely as in [Mo] we use the following : 

Lemma 3.5: For real /U 7^ 0, let = (1 + D/n)~^. Then for large enough, is 
uniformly bounded on L^(R"'"; rdr), and s-lim^^oo Rfi = W. 

Remark: Actually, as in [Mo] we could strengthen the conclusion of the Lemma by replacing 
L^(R+;r(ir) by Hk, k = 0, ±1,±2. Here Hq = L^, 7^2 denotes the domain of a second 
order differential operator (laplacian) on the half-line Aq, = r~^((r9r)^ + rdr — a), a > 0, 
i.e. 7-^2 = {"iA €E : A^i) G L^}, Hi is the closure of Cf^(R"'") for the Sobolev norm 
ll'i/'lli.a = iJ^ drr{\i/)\'^ + \dri^\'^ + :^|V'P))^^^7 and are H-i,H-2 the dual spaces. Actually 
some care is needed because of the singularity at r = 0. We shall not use this however, but it 
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could be also useful, if we want to drop the assumption about the existence of an asymptotic 
of u{r) as r — > oo, to work also on weighted spaces of the form (r)'^L^(R"'"; rdr). 

Given Lemma 3.5, we can prove absence of eigenvalues for An — V as follows. Substitut- 
ing = DR^ for D in the LHS o£(3.18) we get : {(An - V)D^iP\iP) -{D^{An - F)V|V') = 
{D^'ip\{An-V)'ip) = since (An-V) is self-adjoint and D^ip e V{An-V) = L^. Now by the 
definition of i?^, we have [^n — ^5 D^] = — F, D]R^. Using that s-lim^^oo = Id and 

is bounded gives s-lim^^oo -R/uf^n, -Dj-R^u = [^rn-D]. We are left with Rn[AmV]R^. 
Taking r-derivative of (1.10) gives u'{r) = Pf'{PAnu{r)){Anuy{r), and by the explicit rep- 
resentation (1.12) of An 

(3.19) ru'ir) = (27r)-2/3/'(/3A„w(r)) (r[A„, r]u{r) - 2iTnAnu{r)) 

The second term on the RHS is uniformly bounded for r G R"*", since A„ is bounded on L°°. 
We know by Theorem 1.11 that u G C^. According to hypothesis of Theorem 0.2, we shall 
also assume that u has an asymptotics of the form u{r) = nip + ^ + ■ ■ ■ as r ^ cxo. Since 
An has the asymptotic property with vanishing subprincipal symbol, r[A„,r] is bounded on 
m^+X2, so r 1-^ r[Am r]u{r) is bounded on R"*". It follows from (3.19) that the multiplication 
by rV'{r) is bounded on L^, and again by Lemma 3.5, s-lim^^oo Rn[Am V]Rn = [^n, V]. So 
we get {[An - V, Di^]i;\i;) {[An - V, D]i;[ip) as ^ oo, and looking at the RHS of (3.18) 
gives that i/j = 0. 

Absence of singular spectrum follows also from this regularization and Putnam-Kato 
theorem (see [ReSi,Sect.XIII.7]. ) 

By the discussion before Lemma 3.4, the same argument applies to An — W{r), so by 
(3.15), Ln has purely continuous spectrum, which eventually achieves proving Theorem 0.2. 

c) Remarks on the linear stability of higher modes. 

Spectral analysis of operators L^, k > n + 1 is much harder, and we shall content to 
write the formula for Ln+i- For the scalar 1-d Kac model, it is known (see [Pr, Sect. (6.3.1)] 
that if L denotes the corresponding linear operator around an instanton m{x), i.e. a solution 
to Euler-Lagrange equation standing for (1.8) or (1.9), then L is bounded, its spectrum lies 
in R~, and L has a 1-d kernel generated by m'{x). Moreover there is a spectral gap a; > 0, 
namely {'ip[L'ip) < —LuWipW^ for all ip such that {m'['ip) = 0. This relies on the fundamental 
energy identity 

(V'l^V') = -J / dxdyJ{x -y)m'{x)m'{y){ '^\^\ - ^^j^f 
2 J ^m'{x) m'{y) 

which of course doesn't hold in our case. Instead we write, after conjugating with the reflection 

matrix R as before, operator Ln+i in the form : 

~ _ f An+l + An-1 - V{r) An-1 - An+1 \ 

V An-1 -An+1 An+i+An-i-W{r) J 
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The vector U{r) = ("'^f^^)''^) corresponds to the translation mode, and solves Ln+iU{r) = 0. 
With the Gaussian J, naive considerations on the matrix [Ln+i,D] suggest it could enjoy 
positivity properties, but still allowing for a non trivial kernel. 

Appendix 

A) Continuity properties of A^. 

We analyze in this section the action of An on and Y^. 

Lemma a.l: Assume J e *S(R^) (the Schwartz space of rapidly decreasing functions together 
with all derivatives). Let x be a smooth cut-off equal to on [0, 1/2] and to 1 on [1, oo[, then 
if u is continuous on the half-line, we have 

(a.3) sup^<i \ An{l - xhir)\ < Csup^,<i \u{r')\ 

{a A) sup^>ir'=|A„(l < Csup^,<i |w(r')|, k = 0,1,2. 

(a.5) sup^<i \Ar,xu{r)\ < Csup^,>i/2 \u{r')\ 

Moreover, (a.4) holds true for k = 5/2 and for A; = 3 when n > 2. 

Proof. For (a.3), we simply use the fact that Bessel functions (together with all their deriva- 
tives, ) are bounded on R+ (for example Jo{x)\ < 1, and |Jn(3^)| < ^/V2, all n > 1), and 
TJ{p) rapidly decreasing. Next we check (a.4). For A; = 0, the same argument apply. When 
k =1, from the resursion formula for Bessel functions, 

(a.6) (n -h l)Jn+i{x) = x{Jn{x) - J'n+iix)) 

we get 

/■OO /-OO 

rAn{l - x)u{r) = {n+l) dr'r'{l - x)u{r') / dpJn+i{rp)Jn{r' p)TJ{p)+ 

Jo Jo 

nOO nOO 

+ / dr'r'{l-x)u{T') / dpTpJ'n+^{Tp)Jn{T'p)TJ{p) 
Jo Jo 



{a.l) 



By the uniform bound on Bessel functions we just recalled, the first term is estimated by the 
RHS of (a.4). For the second term, we integrate by parts with respect to p in the second 
integral, using that J G 5(R^), which yields the same conclusion. So (a.4) is proved when 
k — 1. Consider at last the case k = 2. Denote by bi{r, r') and h2{r, r') the two p-integrals in 
(a.7). We compute 



(a.8) 



/•OO 7 

rbi{r, r') = {n + l)(n + 2) —Jn+2{rp)Jn{r'p)J^J{p)+ 

Jo P 

(n + 1) / dpJn{r' p)J^J{p) — Jn+2{rp) 
Jo dp 
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To estimate the first term, we spfit the integration over [0,1] and [l,oo], on [0,1] we use 
that Jn{r'p) < Const. (r'p)" since r' is bounded on supp 1 — X- Performing the integration, 
this brings a term 0(r'"). The integral on [l,oo[ can be directly evaluated. We estimate the 
second term in r6i(r, r'), and also rb2{r, r') in the same way. So we proved (a.4) for A; = 0, 1, 2 
and for A; = 3 the same arguments apply when n > 2. 

Let us consider the case k = 5/2. Call 6ii(r, r') and 6i2(r, r') the 2 terms in (a.8), and 
compute 

(a.8) rV26n(r,r') = (n + l)(n + 2) T -%{rpY/^ Jn+2{Tp)Jn{r' p)TJ{p) 

Jo P ' 

Since, for all n 

Jnix) ~ V2{7rx)~^^'^ cos{x — 7rn/2 — 7r/4), a; — > oo 

(rp)^/^ J^_l_2(r/9) = 0{1) uniformly for rp > 0. We split the integral in (a.8) as and . 
Using Jn{f' p) ^ Const. (r'p)"" for r' G suppx, p G [0,1], the -integral contributes for a 
constant times r'^ dpp^~^^'^ = Oi^r'"^), while the Jp^^-integral contributes for a constant. 
Hence r^^'^ dr'r'{l — x)u{r')bii{r,r') < Csup^/^jQ ^^] |'u(r')|. The same arguments show 
this is the case of the corresponding term with r^/^6i2(r, r'), and also of those arising from 
rb2{r, r'). So (a.4) holds true when k = 5/2. 

We proceed analogously for (a.5), and compute A„((-)^'(;)(r), where v{r') = x(r')r'~^w(r'), 
so that |i'(r')| < xi^')^'~'^ ^^Pr'>i/2 I^(^')l7 x(r')r'~"^ is integrable. 

We estimate first ai{r,r') = r' dppJn{rp)Jn{r' p)J^J{p) using (a.6), as r' pJn{r' p) = 
{n + l)Jn+i{r' p) + p-j^Jn+i{r' p) which gives accordingly ai(r, r') = aii(r, r') + ai2(r, r'). In 
0'i2{r,r') we integrate by parts as before, which gives bounded terms for r < 1. Multiply 
again by r', and write a2(r,r') = r'ai(r,r') = a2i(r,r') + a22{r,r'), with a2i(r,r') = {n + 
l){n + 2) ^ J n[rp) J nj^2{f' p)^ J {p)^ and a22{r,r'), a less singular term, containing also 
derivatives of Jn+2{'>^'p)j which we integrate by parts. We multiply a last time by r', use 
again (a.6) and write as{r,r') = r'a2{r,r') = a3i(r, r') +a32(r, r'), with a3i(r, r') = (n + 
l)(?i + 2){n + 3) ^Jn{rp)Jn+sif'p)^Jip)i and 032 a less singular term, containing also 
derivatives of Jn+sir' p). We examine 031 by splitting the integration over [0, 00 [ into [0, r'~^], 
[r'~^,l], and [l,oo]. On [0,r'~"'^], we use the estimates Jn{rp) < Const. (rp)"", and also 
Jn+3{r'p) < Const. (r'p)"^"'""^ when n = 1. So the first integral resulting from that splitting 
is bounded. The integral over [l,oo[ is also clearly bounded. Consider at last the second 
integral, and the most difficult case, i.e. n = 1. This time, we estimate simply J^^^^r'p) 
by a constant, and write Jn{rp) < Const. (rp) ; integrating with respect to p gives a term 
0{rlogr') = O(logr'), but we can still divide this by r'", < a < 1, since r'x{r')r'~^~^^ is 
also integrable. All other terms can be handled similarly. So the Lemma is proved. X 



28 



Now we give some continuity properties of on the space 

Proposition a.2: Assume has the asymptotic property and is positivity preserving. 
Then A^ is a bounded operator on X^, k = 1,2. 

Proof: First wc show that r i-^ Anu{r) is continuous. Let x be a smooth cutoff as above, we 
write Anu{r) = — x)'*^(^) + ^nX'i^(^)- As A„(r, r') is smooth and (1 — x)w of compact 
support, the first term is a continuous function of r. For the second term, write as in the 
proof of Lemma a.l when k = 2, 

(a. 10) AnXu{r)= / dr'r' —f^{r' u{r')) I dpJ='J{p)Jn{rp)r'pJn{r'p) 

Jo f Jo 

and use(a.6) to show that the p-integral is a continuous fimction of (r, r'), uniformly bounded 
with respect to r' e [0, cxo[. Since u G r'^u{r') is also bounded on suppx, and r'^^^ 
integrable, we can conclude that A^u is continuous. 
When k = 1, replace (a. 10) by 

AnXu{r) = J dr'r'—^(r'u{r'))r' J dpTJ{p)Jn{rp)r'pJn{r'p) 

and apply the procedure above. We shall have to estimate Jq°° ^Jn{rp)Jn+2{r'p)J{p), for r 
in a compact set we split again the integration over [0, 00 [into [0,po] and [pojOo]) using the 
bound Jn{rp) < Const. (rp)"^ for the first part, and absolute convergence for the second part. 
Thus we can conclude to continuity as in the case k = 2. 

Next we show that if tt e then r^AnU is bounded near infinity. Write r^Anu{r) = 
r''An{l — x)u{r) + AnXu{r) , by (a.4) we estimate the first term by Const. Hw; For the 
second term, we have for all r > l,|tt(r)| < ||tt; X^||r~'^, so using the fact that An has the 
asymptotic property and is positivity preserving, we get |^ntt(r)| < ||tt; X°||^^(-)~^(r) < 
Const, r^'^lltt; X° II . This brings the proof to an end. X 

Along the same lines, we can prove continuity of A^ in the asymptotic space Y^. 

Proposition a.3: Assume A^ has the asymptotic property and is positivity preserving. 
Then A^ is a bounded operator on Y^, k = 1,2. 

Moreover, if A^ has vanishing subprincipal symbol, then it preserves the closed subspace 
Eo = {ueX^:io{u)=0}. 

Proof: Consider first k = 1. We just need to show that rAnu{r) has a limit as r — > 00 (the 
same as fx, ) and that 

(a.l2) sup,^[,^+^[\r{rAnu{r)-e{u))\<C\\u;Y,^\\ 
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For the first point, we observe that by definition, given £ > 0, — p- < u{r') ^ < p- 

for r' large enough. Multiply this inequality by x(^') before, use that An is positively 
preserving, and An{^)(r) ~ ^, we find [AnX'^i'i^) ~ < v large enough. On the 

other hand, by (a.4) for k = 2, we have An{l — X)'^{''^) = C'(r~^), r — > oo, so we can conclude 
lim^^oo TAnu{r) = £{u). It is also clear that if An has vanishing subprincipal symbol, then 
it preserves the closed subspace Eq = {u E X° : £o{u) = 0}. Now we prove (a.l2). The same 
argument as before shows that if tt e Y^, then \AnXu{r) —£{u)An^{r)\ < C\\u; Y'°||^„^(r), 
so by the asymptotic property, for r large enough, r'^\AnXiJ'{f) — < ^iW + ^ 

C'llwiyiOll. Also by (a.4), we have for r > 1, r^\An{l-x)u{r)\ < Csup^,<i \u{r')\ < C\\u;Y^\\. 
This proves (a. 12) and the continuity property for A; = 1. 

The case k = 2 goes similarly, taking advantage that (a.4) holds for k = 5/2. J(t 

Our next result concerns the behavior of Anu{r) as r ^ 0. 

Proposition a.4: If if -u G X^, then r~"A„(r) is bounded near 0.. 
Proof: We use this time the relation 

2n 

(a. 14) Jn{x) = Jn-l{x) + Jn+l{x) 

X 

p 

which gives, by induction 2^*'(2p)a;~^^ J2p(a;) = ^ ajJ2{p+j){x)i n = 2p, and 2^^'+^(2p + 

j^-p 

p+i 

l)x~^P~^J2p+i{x) = PjJ2(p+j){x), n = 2p + 1. where aj and Pj are rational numbers. 

j=-p 

Then we argue as in Proposition a. 2, for a linear combination of terms of the form 

/•OO i>c>o 

/ dr'r'u{r') / dpp^+'J2p+2j{rp)Jn{r'p)J'J{p) 
Jo Jo 

(we have assumed here n = 2p + 1 the case n = 2p is similar, )and split accoording to the 
partition of unity x ci.nd 1 — X- The latter part is obviously bounded uniformly as r — > 0. 
For the x-integrals, we divide by r'^ as in the proof of Proposition a.2. This easily gives the 
Proposition. ^ 

A more difficult problem would be to analyse the asymptotic behavior of Anu{r), as 
r ^ 0. Formal analogy with radially symmetric solutions of the Ginzburg-Landau equations 
would suggest that if u{r) ~ r"' as r ^ 0, then so does Anu{r). 

B. Renormalized free energy in infinite volume. 

In this section we renormalize (c.31) in infinite volume on the real plane, restricting m to 
a class of functions having topological degree at infinity. As usual, the divergence is due to the 
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self-energy corresponding to a neighborhood of the diagonal in J dxdyJ{x — y)\m{x) — m{y)\'^, 
and to J dx f p{m{x)) . For the latter, we shall remove fp{mp) from fp{m{x)) provided m stays 
sufficiently close to mp. Our renormalized energy needs also to be bounded from below. 

We discuss first some properties of the degree of a complex valued function. See [FoGa] 
for more advanced results. Let m : — > C be a differentiable function, considered as a 
vector field on R^, and subject to the condition |m(a;)| — > as \x\ — > oo uniformly in 
X = x/\x\. Then the integer 

1 /" 1 /" dfTi 
(6.1) degjjm = — / rf(argm) = — / 

independent of R when i? > is large enough, is called the (topological) degree of m at 
infinity, and denoted by deg^^ m. We do not attempt to characterize all functions satisfying 
(b.l), but just state for completeness the following result (the condition on the Vk could 
certainly be weakened. ) We say that m have a one-sided Fourier series if there is n > (the 
case n < follows by complex conjugation, ) such that m{x) = Vn{r)e^'^^ + Vk{r)e^ 



Ake 

k>n+l 



Proposition b.l: Let m have a one-sided Fourier series, Vn{r) — > and m(a;)— v„(r)e* 
as r — > oo uniformly in 9. Assume moreover the Fourier coefficients decay sufficiently fast to 
ensure existence of a continuous limit of Vk{r)z^ as \z\ — > 1~, i.e. k^\vk{r)f' < oo, 

k>n+l k>n+l 

uniformly as r oo. Then deg^^ m — n. 

Proof: By assumption, we can differentiate the series term by term on the circle of center 
0, and radius R. To compute j\x\=R 1^' '^^ P^^ ^ — Since m{x) — fn(r)e*"^ ^ as 

r — > oo uniformly in 9, I z^~'^\ < 1 on |2;| = 1 for R large enough, and Rouche 

k>n+l ^ 

theorem asserts that Vn{R) + ^ Vk{R)z'^~^ has no zeroes inside the unit disc. So the only 

fe>n+l 

pole inside the unit disc is 2; = and the corresponding residue is equal to n. i(t 

Unless restricting to those m with deg^o m = jF(m|m'^) doesn't make sense in the limit 
A —s> 00, so we need a renormalization to remove the logarithmic singularity. More precisely 
we show how to renormalize jF(m), with jF(m) as in (2.1) and m of degree n. For simplicity 
we restrict to the case where m has a single Fourier mode, i.e. m belongs to the sector {en{9)). 
The idea is to remove in (2.1) a thin conic neighborhood of the diagonal x — y m R^ from 
the double integral. 

To give some flavor of the general argument, we forget about the free energy density 
of the mean fleld /^(|m|), and start with the particular case where m{x) = m^e*"^^, which 
belongs to E defined in (2.1) but of course, is not continuous at 0. Denote by T^{m) = 
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I J dxdyJ{x — y)\m{x) — m{y)\^ the interaction term. Let first tq > and spfit 
JpO(m) = l{[ dx [ 

4 J\x\>ro J\y\ 

+ dx dy + dx dy)j(x — y)\m(x)—m(y)\' 

J\x\>rn J\v\<rn J\x\<rn 



(6.9) 



dy 

|>ro 



We are going to reduce (b.9) modulo integrable terms. Because J decays rapidly at infinity, 
it is easy to see that the last 2 terms are bounded. Consider then the first term, and for each 
X in the domain of integration, let L^; = {y e : \y\ > ro, \ty + (1 — t)x\ > ro, Vt e [0, 1]} 
denote the light cone issued from x, and Sx = {\y\ > ro}\Lx the corresponding shadow cone. 
Using the rapid decrease of J, we see that J|a.|>2ro Is d,y\m{x) — m(y)p is bounded. We 
compute I J\x\>ro ■f\y\>ro '^^''^(^ ~ 2/) ^y choosing polar coordinates x = re*^, 
X — y = pe*"^, the volume element is \dx Adx Ady A dy\ = Apr\dr Ad9 Adp Ad(p\, and now we 
perform the integration over 0, = {p > 0,r > ro,(fi e [0,27r],0 e [0,27r]}. Writing y = r'e*^' 
we have the relation 

(6.10) r' sm{e' -ifi)=r sm{e - cp) 

Thus the first term in (b.9) rewrites as 
(6.11) 

- / dx dy J {x — y)\m{x) — rn{y)\'^ — 2rrip I dippdpJ {p)rdrd9 (l — cos n{9 — 9')^ 

J\x\>ro -'|a;|>ro JO, 

We make the following observations : for fixed let be the reflection on the line 9 = (f. 
Then for all r > p, the map {r,9) ^ {r',9'), is 27r-periodic in 9 (it corresponds to shifting 
the circle of center and radius r by the vector p(cos(/?, — simp). ) For given r, the function 
^ I— r' is even under F,^, and increases from r — p for = to r + p for = (/? + tt, while 
9 ^ 9' — 9 \s odd under F,^, increases from for 9 = p to a, maximum value {9' — 0)max, with 
cos((0' — 0)max) = ^(r^ + /9^)~^/^, lox 9 ~ p) + 7r/2, and decreases again to for 6* = <y? + tt. 
Together with relation (b.lO) this gives 

(6.12) \r' - r\/r < p/r, \{9' - 9)^^^\ ^ p/r 

We compute (b.ll) as follows. For C > large enough to be fixed later, we split the integral 
over the sets {r > Cp}, and {r < Cp}. For the first part, {9' — ^)max is a (non degenerate) 
critical point, and we observe that for C > large enough : 

(6.13) 2 / d9{l - cosn{9 - 9')) = 8 ^ d9{l - cosn{9 - 9')) = TTn^{^f + 0{^f 

Jo Jip J. J. 
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3 3 

where is asymptotic, as p/r — > 0, to Cn(^) , for some Cn > 0, all n (cn increases to 

+00 with n. ) This shows that the corresponding integral contributes with a logarithmic 
singularity 



/ 



(b.U) 7rn^ / dpp^J{p)y, no = {r>ro,r>Cp} 



which we substract from 2 J^d9d(prdrpdpJ{p)(l — cosn{d — 6')). [Actually we replace the 
integration over r > ro in (b.l3) and (b.l4) hy N > r > ro and let then N — > oo.] Since the 
^-integral is clearly independent of (p, integrating over (p gives an additional factor of 27r. The 
term C(^) in (b.l3) then contributes to a finite, and positive integral. Consider next the 
integral over r < Cp. Interchanging the drdp integrals (which is legitimate after the cutoff 
r < N, ) we are lead to estimate 



/•oo /"OO /•2tt /•2tt 

2 rdr dppJ{p) dcp de{l - cosn{e - 6')) 

Jo Jr/C Jo Jo 



There we use simply that the dcpdO integral is bounded and take advantage of the rapid 
decrease of J to bound the p-integral by a negative power of r to make convergent the 
resulting r-integral. Thus we proved 

Lemma b.4: With the notations above, if m{x) = m^e*"^, let 

ren(- J dx J dyJ{x — y)\m{x) — m{y)\^^ 

= - J dx J dyJ{x — y)\m{x) — m{y)\'^ — 2{'ri:nmi3)'^ J dpp^J{p) — 



(6.15) 



Then < < ren(|; J dx J dyJ{x — y)\m{x) — m(y) < +oo. 

Now we turn to the more general case, which is sufficient for our purposes, where m{x) 
belongs to the (eTj,)-sector, and |m(x)| mp as x — oo. First we renormalize the mean field 
free energy jF^(m) = J dxfpimix)) for \m\ — rrifj e L^(R^) ; taking advantage of the fact 
that //3 attains its minimum at m/3, we let simply 

(6.16) < J^reJm) = J dx{ff3{m{x)) - fp{mp)) < +00 

he first inequality because . Next we pass to the interaction term J-'^{m) as we did for 
m{x) = mpe^'^^ . So let ro > and split the integral as in (b.9) ; as before we are left with 
the integral over {\x\ > tq, |y| > ro}. Assume also that m is continuously differentiable 
there. Use Taylor formula to rewrite m{x) — m{y) = Vm{ty + (1 — t)x)dt ■ {x — y), let 
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z = ty+ {l — t)x = re'^^ ,x — y = pe^"^ . Expanding the product, multiplying by J{x — y) = J{p), 
computing the Jacobian, and integrating, we find 

^ r r poo /•27r 

- / dx dyJ{x — y)\m(x) — m(y)\^= dt dt' / J{p)p^dp I dtp 

4 i|a;|>ro J JO Jo Jo Jo 

dej^ rdr[—{z) — {z')cos{<p-d)cos{<f-e') 



(6.20) ^ 



TO 

+ ^ ^ (^) S (^') «-(^ - ^) - + ^ (-) i S sin(^ - e') cos(^ - 6) 



r de^ ' dr'^ ' ' ' dv^ 'v' dd' 

1 dm 1 dm 



where z' = t'y + (1 ~t')x = r'e^^ . Because y G L^;, the segment [x, y] lies outside i?2(0, tq) so 
r, r' > ro and all terms in the integral on the RHS of (b.20) are well defined. To start with, we 
make as before a cut-off in the r-integral in the RHS, so that we can split it into different terms. 
The domain of integration is denoted by f2 = {0 < < 1, r > ro, p > 0, (/?, 6' G [0, 27r]}. 
The observations leading to (b.l2) can be exactly repeated, changing x to z, y to z' , and p 
to \t — t'\p, now we get : 

(6.21) \r'-r\/r < \t-t'\p/r, cos((^'-^)„,ax) = r(r2+(t-t')'p')~'/', |(^ '-^)max| ~ \t-t'\p/r 

Making use of the symmetry in 9, 0' , (b.20) is the sum of 3 integrals, whose corresponding 
integrands write : 

A{t, t\ p, r, Lpfi) = u'{r)u{r') cos(99 — 9) cos((/? — 9') cos n{9 — 9') 

B{t, t\ p, r, ^p,9) = 2n[^^u'{r') sm{(p - 9) cos{^p - 9') 

(6.23) uir') 

+ u'(r)^;^ sin(v; - 9') cos{^ - 9)] smn{9 - 9') 

D{t, t', p, r, ip,9) = V? "'-^-'y^' sin((^ - 9) sin((^ - 9') cos n{9 - 9') 

rr 

{u'{r) denotes the r-derivative of u. ) We examine first D, which we renormalize by setting 
u{r) = mp + v{r) ; using (b.lO) we find : 



(6.24) D = [(^f + "^Hr) + v{r')) + r-\{r)v{r')\ sirv^i^ - 9) cosn{9 - 9') 



Comparing the coefficients of m| in (b.20) and (b.ll), (b.l3) shows that we must have, for 
r > Cp : 

(6.25) ""^^r^^ Iq - ^) - = ^^^(^)' + ^{^f 
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Since the ^-integral is clearly independent of integrating over ip gives again an additional 
factor of 27r, so that J~dtdt'd(pJ{p)p^dprdrn^(^^Ysia^{ip — 9)cosn{9 — $') contributes to 
(b.20) [after letting A?" — > oo] with the logarithmic singularity 

2{7rnmpf dpp^J{p) — , = {r > tq^t > C p} 

which we eventually substract from (b.20), and the term C(^) in (b.25) then contributes 
with a finite, and positive integral. 

Consider next the integral over tq < r < C p. Interchanging the drdp integrals we are 
lead to estimate 



/■oo fco i'2iT /•2'ir 

/ drr{^Y / dpp^J{p) / dip ddsm^{ip-d)cosn{e-d') 

Jro Jr/C Jo Jo 



There again, we take advantage of the rapid decrease of J to bound the p-integral by a 
negative power of r to make convergent the resulting r-integral. Now we examine all other 
terms contributing to (b.20), namely A{t, t', p, r, cp, 9), B{t, t', p, r, (p, $), and the part 

D'(t, t', p, r, ip, 6) = ri" [^{v{r) + v{r')) + r-\{r)v{r')\ siTL^{ip - 6) cos n{e - 6') 

which is left from D. Consider first D', and recall v e L^(R"'"; rdr). For the first term we use 
Cauchy-Schwarz inequality to get 



l*OG POD POD 

/ rdr^^;(r) sin2(^ -^) cosn(^ < (/ rdr\v{r)fy^{ rdr{^fy^ 



and similarly for the second term, if we think of the fact that r and r' play symmetric roles. 
Integrating against J{p)p^ with respect to dtdt'd(pd9dp gives finite quantities. The last term 
in D' contains the correlations r~'^v{r)v{r'). Again we split the integration according to 
{r > Cp} and {r < Cp}, and for the second part, use the rapid decrease of J. For the 
first part, we use instead that the "translations" (t, t', (/?, ^, p) i— > (r i— > f(r')) are uniformly 
continuous in L^([Cp, +oo[;rdr) when Cp < tq ; this gives 

I / dtdt'dipJ{p)p^dprdrd9r~^v{r)v{r') 

jQ,r>Cp 

X sm'^{(p-e)cosn{d -e')\ < Const. \\v; L'^{[ro,+oo[;rdr)f 

Consider then A or B. These terms involve the derivatives r~^v'{r'), r'~^v'{r), and the cor- 
relations v'{r)v'{r'), r~^u{r)v' {r'), r'~^u{r')v'{r), so they can be treated as above, provided 
r^v'{r) is in L'^irdr) for some 5 > 0. This holds true when u{r) ~ mp + C(l/r), r — > oo and 
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this relation can be differentiated, e.g. if m has the asymptotic properties of a symbol in 1/r. 
Summing up, we proved the : 

Proposition b.5: If m{x) — u{r)e^'^^ is bounded, u — + f , with 

V e W = {L^([0,l];r"^dr)niy^(R+;rdr),(-)V e ( [1, +00 [; rdr)} 
for some 5 > 0, then : 

■^ren.roM =\ J J dyJ{x - y)\m{x) - m{y)\'^ - 2{TTnmpf j dpp^J{p)y 
+ J dx{f(3{m{x)) - fp{m(3)) < 00 



(6.28) 



Moreover, J^^en is (strongly) continuous on W endowed with its natural Hilbert space struc- 
ture. 

Remark b.6: The renormalization can be easily extended to the case where u has a disconti- 
nuity on a sphere r = A, i.e. when it has derivatives a.e. This is used when considering the 
partial dynamics. 

When m is bounded, it is clear that .7^°en,ro ("^) -^njn.n ("^) '^i^^ J^^^ differ by a 
quantity of order (tq — ri)^. Next we show that if J > 0, then .Fren,ro("^) is bounded 
from below in some region of the configuration space. First we consider the interaction 
term .F°(m) = \ J dxdyJ{x — y)\m{x) — m{y)\^, and compute the formal Fourier transform 
of /i : R2 ^ C, h(x,y) = y/J(x-y){m(x) - m{y)). We have J^h(^,r]) = {j^Vj(-r]) - 
J^Vj{^))j^m{^ + r]), so by Parseval identity (still in the formal sense) 

I dxdyJ{x-y)\m{x)-m{y)f = {27r)-^ J d^dr]\J^m{^ - ri)\^ {j^y/j{r]) - Ty/Ji^f 

so we have exchanged the roles of m and J. Of course the integral is divergent because 
m ^ L^, but the singularity is due to \ J-'m{^ ^^)Pi ^lot to the boundary condition at infinity, 
since deg^ J^\/j = 0. The finite part (P.f.) of |JFm(^ — 77)^ is |J^m(^ — ?7)p|g=r, = (/ mxlx)'^ 
and is finite when m{x) = e*"^w(r), u{r) = rriis + v{r), v e L^, because of the periodicity in 
9. So another renormalization of J^^{m) is given by 

•^rWM = ^(27r)-4 Pf. j didr^\rm{i - r^)\^{r^{n) - TVJ{0)' 

which is obviously positive. Moreover, Lemma b.4 and its proof show that both renormaliza- 
tions should agree up to a finite term, depending only on tq, when m{x) = mpe^'^^ . Thus, we 
should have .?^ren,ro ("^) = •^ren,:F(™) niodulo a constant term, depending on ro, but not on 
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m, when m satisfies all hypotheses of Proposition b.5. In particular, ^°en,ro("^) bounded 
from below. 

C). The free energy in Kac's model with continuous symmetry. 

In this Section, we recall and make more precise the procedure of renormalization in the 
continuum, carried out in [AlBeCaPr], [Pr] for spins valued in {—1, +1}, and extended to the 
XY model in [BuPi], when a vorticity condition holds at infinity and the interaction is not 
necessarily compactly supported. Actually, the final form for the continuous renormalized 
free energy functional should be regarded as a postulate. 

Note that passing from the lattice to the continuum amounts to consider convergence of 
Riemann sums as the mesh of the (scaled) lattice goes to 0, as an homogenization process, 
i.e. the convergence of discrete measures. 

a) Some definitions. 

Consider the lattice Z'^, consisting in a bounded, connected domain A (the interior 
region), and its complement (the exterior region) A^ (twiddled letters will always denote 
discrete objects on Z'^. ) Physical systems make sense in the thermodynamical limit A Z'^, 
in the sense of Fisher. The simplest way of taking this limit is to double the side of the unit 
hypercube repeatedly, so the side of A = A~is of the form 2^, £ G N. 

To each site z G Z^ is attached a classical spin variable a{i) G S"^. The configuration space 
A'(Z^) = (S'?~^)'^ is the set of all such classical states of spin ; it has the natural internal 
symmetry group 0^{q) {q = 2 for the planar rotator.) Given the partition Z'^ = A U A'^, we 
define by restriction the interior and exterior configuration spaces A* (A) and A'(A'^), and the 
restricted configurations by a~ and cr-^. 

It is convenient to rescale Ai to a domain of fixed size L = 2^ £ G N, A c R^, with 
scaling factor is 2~^ ; for the moment we could think of A as the unit square {£ = 0) but we 
shall eventually let also A ^ oo. 

Following [Pr], for /c G N, we denote by Q^'^^ the partition of R'^ into small cubes 
(j{k) ^ 1^ ^ (ri,---,rd) G 'R'^,2~'^Xi < Vi < 2~^{xi + 1)}, of side 2~^ , and indexed by 
X = {xi, - ■ ■ ,Xd) G Z*^, called (Q*^'^)-)atoms. The atom C^'^\r) is the unique atom of Q^^^ 
that contains r. We say also that a function on R'^ is Q'^^-'-measurable if it is constant on 
each atom of Q^^\ and a set A C R'^ is Q'^'^^-measurable if its indicator function is Q^^^- 
measurable. This allows in a natural way to identify a function a on the lattice with a 
Q(fc)_nieasurable function a^'^^ on R'^, assuming (j^^\r) — cr(x) with x — {xi, ■ ■ ■ , xa) and 
r = (ri, ■ ■ ■ , rd) as above. 

Let now 7 of the form 7 = 2"^^^ , G N, that will be the inverse of the interaction length 
in Kac's potential. Given a state a G we define as the Q*^'^'>')-measurable function 
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crC^T). Since we take a simultaneous limit A — > oo, 7 — > 0, rather than Lebowitz-Penrose limit 
A ^ 00 followed by 7 ^ 0, it may be convenient to label the configurations by 7, instead of 
A. We call also cr^ a smooth- grained configuration. 

Because of the internal continuous symmetry of A" (A), the probability distribution v 
for the states of spin is defined as the normalized surface measure on S^~^, i.e. ^{dai) = 
Lij~^5{\ai\ — l)dai, where ujq is the volume of S^~^. 

Given as above, and an integer Uj < kj, we associate the Q^"'"') -measurable function 
(magnetization) 

(c.l) m,{r) = 7r(--)^^(r) = / dr'a,{r') 

These averages of a-y over the "intermediate" boxes (or block spins) C^'^^\r) of volume 
7'^'^, define the coarse-grained configurations, and the map tt^'^t) : 5^(0,1), = 

(S«-i) , 5^(0, 1) the closed unit ball of R«, is called the block -spin transformation. We 
may think of it as a random variable. 

N 

More generally, let ttn : (S*?"^) B'g{0, 1), a = (cti, • • • , ctat) ^ 7:N{cr) ^ jf"^ (Ji- It 

is easy to see that ttn is a smooth, surjective map, and its restriction ttn '■ (S^"-*^) \ A — > 
Bq{0,l) a submersion, where A denotes the diagonal (7i = • • • = o-jv of (S«-i) , andSq(0,l) 
the open unit ball of R^. Hence the probability measure i/n = i'^N)*ii^ <S> • ■ ■ <^ 1^) has 
a smooth density with respect to the (normalized) Lebesgue measure on Sq(0, 1), namely 
^^(m) = J/ _ N JV YliLi i^{da{i))5{7rN{a) — m), and we can check as in [BuPi,formula (5.1)] 

that 

dl^N , X /N.q f , -iNlv.m.\/ f /, \ ih,.fr\\N 



dm 



(here (■, ■) is the standard scalar product in R^. ) 

Of course, Kolmogorov realization theorem would allow to define this way a probability 
measure on (-Bg(0, 1))^, but this will not be used at the present level. 

Now, following [Pr] wc choose rij so that 7'' = 2""^ for some < 5 < 1, a good choice 
is 5 = 1/2. From the point of view of the discrete scheme, a coarse configuration is defined 
on a "coarse lattice" A* which we magnify by the factor 2^~>~'^'^ to the "smooth lattice" A 
(see e.g. [El-BoRo]. ) Thermo dynamical properties of the system are most significant on the 
coarse lattice. 

From the discussion above, for /c G N, let C^''^ = 2^C^^'> n Z"' denote the atom C''^^ 
rescaled to the lattice units ; thus, C^^'^ is the rescaled block-spin. So the probability distri- 
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bution of the empirical average tt^'^t) has Radon- Nikodym density 



We shall sometimes write (T-y{i) instead of cr~(z). Let also N — \C^'^'^^ \ = ^("^"i)*^, ;/("-t) = ^y^. 

The continuous Kac Hamiltonian is defined as follows. Let J > be the interaction 
potential, for a given a'^ e X{A'^), and for A = [jC^'^'f^ = [jC^'^'^^ as above, we define the 
energy of a configuration a e '^(A) by 
(c.3) 

H{a^\a;) = -l [ dr [ dr'J{r-r'){a^ir),a^{r'))- I dr I dr' Jir - r'){a^{r), a'^{r')) 

^ J A J A J A JA'^ 

Again from the point of view of the discrete scheme, since c^^(7^) = cr(z), for i e Z'', there is 
a corresponding Hamiltonian on defined by 

(c.4) H^{a~\a~^) = -]-^ J^{i,j){a{i),a{3))- ^ J^{i,j){aii),a{j)) 



2 

i,j€A (iJ)eAxA= 



where 



(c.5) J-yiiJ)^!'"^ [ dr I dr'J{r-r') 

The two hamiltonians are simply related by : 

(c.6) #7(a~|a~ )=7-'iyK|cT^) 

so that H is an intensive hamiltonian. As observed in [AlBeCaPr] , neglecting the variations 
of J in the integral, we get 

(c.7) Mhj)^l^Jil\i-j\) 

which has the typical scaling properties of the original Kac potential, and the results of this 
Section remain valid when the energy is given by (c.4) with (c.7) holding as an equality. 
This observation allows to define Gibbs measure conditioned by a~^, on the space of spin 
configurations X{A) with mesh 7, at inverse temperature as 

(^•S) /^/3,7,a('^^aI^ac) = z — exp[-/3^7((7~|c7~ )] II u{da^{i)) 

/3,7,A^%^^ ieA 
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where 

(c.9) Z^^^~{a~J= [ l[i.{da~ii))exp[-pH^{a~\a~J], Qo = {S'^-'f 



is the partition function, making of 1^ p ^ 'j^i'^^W^c) ^ probabihty measure on the product 
space. As before, we can take the direct image //^ ^ a^'^'^aI'^A'^-^ through the block-spin trans- 
formation. We shall discuss this in the following subsections. 

b) Entropy estimates. 

We want to relate jj-^og ^^(m) with the entropy functional I{m) of the mean field 
approximation. Recall I{m) = sup^^gj^, ((/c,m) -log0(/c)), (^(/c) = /s,_i e<'='">zy(d'y), and 
/(m) = {k*,m) — log(p(k*) where k* — k*{m) is the unique point in R"^ that achieves the 
maximum. Clearly also, by spherical symmetry, /(m) = /(|m|), and = sup^>Q(t|m| — 

log^(t)) ^ t*\m\ - log0(r), where t* = t*{\m\), t*{0) = 0,r(p) ~ (2 - 2p)-^ as p ^ 1. 
Furthermore the supremum is achieved when k* and m are colinear. For |m| < 1 we introduce 
the probability measure iJ,{dx;m) on S^~^ defined by 

(c.ll) iJ,{dx;m) —exp{{k*,x) — log(f){k*))i'{dx) 

As in (c.2) we define 



dm 



J . — - 



and denote HN{dx;m) = p^'^'^\dx]m). It is easy to see that ^^("^) = e ^^^^^ ^^^^{'m)- Let 
also be the complex function defined on R"^ 

Recall from [BuPi,formula(5.7)] the identity 

The observation is that -^^og^^{m) is a small correction to —I{m) (the entropy for the 
mean field) as N becomes large. Indeed we have : 

Lemma c.l: Let g = 2 for simplicity. With the notations above, we have for all A'' > 1, 
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for some Cq > 0,q' > 0. 

Proof. We need to show that the integral on the RHS of (c.l2) grows at most hnearly in 
(1 — |m|)~^/^, with coefficients polynomial in N. Notice first that < v'm(O) = 1 for 

all V e R^, so integrating ipm{fj)^ over the ball in of center and radius N we get 

(c.l5) \[ dv(pm{v)^\ < Const. N'i 

J\v\<N 

Now, we estimate the integral near oo, using complex stationary phase. We will be a little 
sketchy, but it is easy to see that our leading terms give the correct behavior with the required 
uniformities (see e.g. [Sj] for more details. ) Let v = r(cos<^, sin(^), (p e [— tt, tt], we rewrite 
<Pm{v) as 

(c.l6) (^^(^;) = (27r0(r))"^e*^l^l^'°'^ / d9e-'''^^^''^\ ${t*) = loit*) 

J —n 

with ^{9,(p) = cos(^ — (p) + zAsin^, and A = t*/r. Here we consider r > N as the large 
parameter. The critical points in 9 are given by the equation sin(^ — (p) — iXcos9 = 0, 
so ^ I— > ^{9, (f) has no real critical point if ^ iTi'/S, and 2 real critical points 9 = ±<^ 
otherwise. Actually, Im$(±7r/2, ^7r/2) < so the contribution of the critical point with 
sign opposite to this of will be exponentially small, and by symmetry it suffices to consider 
{9,ip) = (7r/2,7r/2). This is a non degenerate critical point, since ^^(7r/2, 7r/2) = — 1 — iX. 
Because of analyticity, there is a complex critical point 9c = 9c{<p) for nearby values of (p, 
and a simple calculation yields 

(c.17) $(^e, ^) = 1 + iA - ^^YTo) " ""Z^^' + " ""Z^^' 

where 0{<p — 7^/2)^ is uniform in A. The complex Morse lemma then shows that the local 
analytic diffeomorphism 6 ^ 6 given hj9 = h{9 - 9c;<p)e-'^/'^{9 -9c), /i(0; 7r/2) = 1 + iA, 
is such that $(0, = $(^C5</') ~ 19"^ / 2. Then complex stationary phase shows that the 
contribution of a (fixed) neighborhood of 9c to the integral in (c.l6) is given at leading order, 

1 /2 

by ( yj^^(o.y,) ) e*'^/^e~*'^*''^= Outside this neighborhood, non stationary phase arguments 
show that the integral is exponentially smaller than e~'^^'^'^^'''^\ as a function of r, and we 
eventually get, for in a (real) neighborhood V± of ±7r/2 : 

(c.l8) r d9e-'^^^^'^^ = ( j '^V"/V^"'^^^-^)(l + i?i(y^,A,r)/r) 

with \Ri{ip, A,r)| < Const.. By analyticity, we can keep track of the critical point 9c for all 
ip e [— 7r,7r], so formula (c.l8) still makes sense when p) ^ V±, but then Im$(^c)¥') < ■, 
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and ^ V± doesn't contribute to the final result. Now we raise (c.l6) to the power N and 
integrate over r > N, using (c.l8) this yields (with a factor 2, accounting for the contribution 



[ dv<pm{vf ^ 2(27rJ(r))-^e^^'^/'^(27r)^/2 / rdr{rh{0; ^))-'"^ 

^|z;|>ri Jr>N 



where $i(</7) = ^{9c,f) — \m\ simp. Then (c.l7) shows that $i has a non degenerate point 
at (f — 7r/2, and as before, the complex Morse lemma gives a local analytic diffeomorphism 
cp^ip = f^{cp- 7r/2)(v? - 7r/2) with §{n/2) = ^(O) = + i\m\Y^\ such that $i(v?) - 

$i(7r/2) — i^^, $i(7r/2) = 1 — \m\ + iX. So by complex stationary phase, with Nr as a large 
parameter, we can evaluate the inner integral in (c.l9) : 
(c.20) 

Jv+ 

^_,jvr(i-|m|)givr y^^^/2(0)-^(l + i?i(7r/2. A, r)/r)^(l + i?2(A, |m|, r)/iVr) 

with i?2(A, |m|,r) < Const.. At last, we estimate the resulting r-integral in (c.l9). This time 
we compute simply an upper bound for the integrand. We have (1 + Ri{ti /2^\,r)/r)^ < 
Const, uniformly in as r > A'^, and the same holds for (1 + i?2(A, |m|, r)/Nr), so inserting 
(c.20) into (1.19), we get 

(c.21) / dvifmiv)^ < Const. [ rdr{^J^—-\r + it*\-^/^)^^/2^^\f2{0)\-^ 

J\v\>ri Jr>N V27r/o(t*) 

We have 



and since ,^ f ^ , , ^ 1~ as t* — > +oo, it is easy to see that 

v27ri*Jo(r j 

limsup(^J^ (A2 + r2)-i/4)^ < Const. 



uniformly in t* > 0. Furthermore |/2(0)|~^ = (r^ +r^)^/^(r^(l - \m\)'^ +r2|mp)~^^^, and 
since t* ~ 2~^(1 — Iml)"-*- as |m| — >• 1, we have |/2(0)|~^ < Const. (1 — \m\)~'^/'^ uniformly as 
r > N. So the integral on the RHS of (c.21) is bounded by a constant times 

A2+r2 



(I-H)-VW rdrV2^rC^^±^) 

Jr>N ' "T 
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and it is easy to show that there is q-' > such that j^^j^rdr^j2TT/Nr{ ^2^l*i )^^^ < 
Const. N"^ uniformly in t* > 0. This completes the proof of the Lemmc. X 

Next we have the following : 
Proposition c.2: With the notations above, there is Ci > such that 

(c.25) I log + 7-'^/(m, A)| < C^{L^-^f\og^-^ + log J] (l - |m(a;)|)"'/' 

a;eA/C("-T) 

Proof. Sum (c.l3) over all the cubes C^^^^ contained in A, which have same cardinal N = 
^{6-i)d^ and multiply by N the resulting equality. The first term on the LHS 



(c.26) log n -^{m)^ ^i.{d<j~(i))5{ir^-'^^<j~(i)-m) 



is the Radon- Nikodym density ^^^(^) of a probability distribution, namely the family of 
empirical averages tt*^""'') in C^'^i^ considered as i.i.d. random variables. Here we need interpret 
m as a Q'^'^T^-measurable function on R'', and the RHS of (c.24) should actually read 

where somewhat incorrectly, the notation A/C^'^'<^ reminds us we have tiled A by translates 
of C^'^^\ and ci""^^ is the atom of the partition of A labelled by x. We have IA/C^^t^] = 
|A/C("'t)| = {L^~^Y. This identification will be made freely in the sequel. Summing up 
the second terms on the LHS of (c.13) over x will produce 7~'^/(m, A), where /(m. A) = 
drl{m{r)) defines a functional on the space of Q^'^T^-measurable functions. While summing 
up the RHS of (c.l3) over a;, (c.25) follows easily if we make use of the inequality log(a + 6) < 
log(l + a) + log 6 valid for any a > 0, 6 > 1. X 

c) Free energy estimates. 

Following [Pr], wc replace now in the hamiltonian the spins C7(r) by the magnetization 
m-y(r) as in (c.l). We have : 

Lemma c.3: With the notations above, for some C2 > we have : 
(c.30)) \H{m^\m't^) - H{a^\a''^)\ < C2LVIIVJII1 

(here || • ||i denotes the L^-norm on R*^) 
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Proof. By the definition of H(aj\a^) and (c.l) we have, denoting C^(r) for C^'^-'^r) 

H{my\m^) - H{a^\a^) = (l- [ dri [ dr2 + [ dri [ dr2){(T^{r\),a^{r2)) 

^ J A J A J A Ja= 

X [7-2'^'^ / dr [ dr'{J{r-r')-Jiri-r2))] 

where we have used ri e C,y{r) iff r e C^(ri). We estimate J(r — r') — J{ri — r2) by Taylor 
formula, noticing that |(r — r') — (ri — ^2)! < Const. 7*^, so the dr2 integrals over A or A'^ 
are bounded by a constant times 7'^||VJ||i, then the resulting dri integral over A gives an 
additional factor, which proves the Lemma. X- 

Introduce the free energy in A at inverse temperature /?, inclusive of the interaction on 
A*^, as the functional on the space of Q*^"'') -measurable functions m(r) 



(c.31) 



T{m\rrf) = \ [ dr [ dr' J{r - r')\m{r) - m{r')\'^ 
4 Ja ,/a 

+ 1 [ dr [ dr'J{r-r')\m{r)-m{r')\''+ [ dr{fp{m{r))-fp{mp)) 

^ J A Ja<= J a 



where we recall fpim) = — ||mp + ^I{m). Let e be any (fixed) unit vector in R^, and 
the constant function on A equal to m/je, which we extend to be equal to m*^ on A'^. The 
functionals /(^A) and H{-\m^) are related to J^{-\m^) in a simple way : 

(c.32) .F(m|m^) - J^im^lm") = {H{m\m'') + ^/(m. A)) - {H{rhi3\m'') + ^/(m;3. A)) 

Analogously to (c.8), (c.2) we introduce the canonical Gibbs measure conditioned by the 

external configuration = : 

(c.33) 

l^p,y,A{day;m\a^) = -— / exp[-P^-'^Hy{(7^\a^)] i^{da~{i))6{7r^(Ty{i) - m) 

^0,1, AW-) Jn„ _ 



i£A 



where the partition function Zf^^^^Aia^) = ^/s^^^'^A"^ defined in (c.9). We have made 
use of (c.5) to work on the rescaled lattice A, and set 5(7r-),a"-y(z) — m) = 5(7T^'^''^a~{i) — rti). 
By definition of the image of Gibbs measure through the block-spin transformation, we have 

(c.34) / (im//;3,-y,A((ia"-y; m|a"^) = 1 

J\m\<l 

where dm is the normalized Lebesgue measure on the product space HxeA* ^q(^^ Next 
we give a precise meaning to the approximation A*/3,7,a ~ exp[— /37~^^(m|m'^)] stated in the 
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Introduction, by establishing the analogue of [AlBeCaPr, Lemma 3.2] in case of continuous 
symmetry, improving also [BuPi, Lemma 3.1] : 

Theorem c.4: With the notations above, there are constants Ci,C2 > such that for any 
coarse-grained configuration m on A, we have : 

-(7(?Ti) - (L7-i)^(C2/37iVJ||i + Ci7(i-'^)'^log7-i) 

(c.35) < log iJ,i3^j^A{daj; m\a^) + l3^~'^T{m\rrf') 

< g{m) + inf^^g, T{mp\rrf) + {L-t-^)\C2(5i^\\V J\\i + C^-^^^-^^'^log-f-^) 

where g{m) = log n^eA/c{"7) ~ \'m{x)\) 

Proof. First we look for a lower bound on Z^^^^a(c^)- Using (c.30) we get : 

exp[-/37-'^iy((7^|a;)] > e^^[-p-f-'^H{TT^a^\a''^)] exp[-C2/?7"''^Vll V J||i] 

(and similarly for the upper bound). Multiply these relations by 5(7r^(j^(i) — m), integrate 
with respect to Jlj^A ^('^'^A^^-^) ^^^^ (c.25), we get 

exp[-/37~^(iy(m|m'') + -^/(m. A))] exp-7/;^(m) 
(c.36) ^ / n^(^'^A(^))®^P[~^^"'^^KI^7)]^K^7W -"^) 

i€A 

< exp[-/?7~'^(iy(m|m'') + -|/(m. A))] expV'7(m) 

where 

(c.37) = (i^7~')''(C2/37'l|VJ||i + Ci7^'-'^'^log7-') +^(m) 

Now we estimate the contribution of a neighborhood of rhp to the partition function. So let 
< X < 1 be a smooth positive cutoff equal to 1 near m = m^j, multiply (c.36) by x(™) ^'^^ 
integrate over m with respect to the product measure YY^e'^/c'-^-y^ dm{x), then make use of 
(c.9), (c.32) and (c.34), we get 

X y (imx(m) exp[— /37~'^(jF(m|m'^) — J-'{mj3\m'^)] 

where HV'ilIx ~ ^^Pmesuppx ^7("^) ^ Choose supp x so small that \m — < 
on supp X- Using (c.31), the normalisation of J and Taylor expansion of around m/3, we 

on supp X) so : 

(c.38) 

^/3,7,a(^^) > exp[-(||V^^||^ + C3(L7-i)^/37(i-^)'^)] exp[-/37-'(Jf(m^|m'^) + 
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Inserting this in (c.33) and (c.36) we find 
(c.39) 

For the upper bound on Zp^y^\{a^), we use (c.34) and (c.36) to write 

J J ^ 

< J dmexp[-P^~'^{H{m\m'^) + ^/(m, A))] expV'7(m) 
By (c.32) and inequahty J^{m\m^) > 0, we have 



X J dm exp ijj^ (m) 



se we are left to estimate J dm exp i/)^ (m) , the integral running over the product space, and 
'0-y(m) as in (c.37). Since /q^(1 — p)~^^^pdp < oo we find 

^/3,7,a(^^) < exp[-/37-'^(i/(m^|m'^) + ^/(m^, A) - .F(m^|m^))] 

X exp[(L7-i)'^(C2/?7iVJ||i + Ci^^^-^^'^log^-'' + C^i^^'^^'^)] 

Inserting this and the first inequality (c.36) in (c.33), we find, absorbing the C4-remainder 
term into the Ci -remainder term, and using 

\ognp^^^AidaY,m\a^) > -g{m) - /3-r'^\F{m\m'') 

(c.43) 

- (L7-^)^(2C2/37'^||VJ||i + 2(717^1-'^)'^ log 7-1) 

Putting (c.43) together with (c.39) with new constants Ci,C2 gives the Theorem. ^ 

Of course, these estimates break down when |m(x)| gets close to 1 for some x E A, which 
reflects the fact that the entropy density I{m) is singular near \m\ = 1. It is shown in 
[BuPi, Theorem 2.2] that z^Ar({|TO| > 1 — p}) decays exponentially fast as — > 00 when p > 
is small enough. 
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